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ABSTRACT 



When a quantum field theory has a symmetry, global or local like in gauge 
theories, in the tree or classical approximation formal manipulations lead to believe 
that the symmetry can also be implemented in the full quantum theory, provided 
^ . one uses the proper quantization rules. While this is often true, it is not a general 

^ ' property and therefore requires a proof because simple formal manipulations ignore 

the unavoidable divergences of perturbation theory. 

The existence of invariant regularizations allows solving the problem in most 
cases but the combination of gauge symmetry and chiral fermions leads to subtle 
issues. Depending on the specific group and field content, anomalies are found: 
obstructions to the quantization of chiral gauge symmetries. Because anomalies 
take the form of local polynomials in the fields, are linked to local group trans- 
formations, but vanish for global (rigid) transformations one discovers that they 
have a topological character. 

In these notes we review various perturbative and non-perturbative regulariza- 
tion techniques, and show that they leave room for possible anomalies when both 
gauge fields and chiral fermions are present. We determine the form of anomalies 
in simple examples. We relate anomalies to the index of the Dirac operator in a 
gauge background. We exhibit gauge instantons that contribute to the anomaly 
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in the example of the CP{N — 1) models and SU{2) gauge theories. We briefly 
mentioned a few physical consequences. 

For many years the problem of anomalies had been discussed only within the 
framework of perturbation theory. New non-perturbative solutions based on lattice 
regularization have recently been proposed. We describe the so-called overlap and 
domain wall fermion formulations. 
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1 Symmetries, regularization, anomalies 

Divergences. Symmetries of the classical lagrangian or tree approximation do not 
always translate into symmetries of the corresponding complete quantum theory. 
Indeed quantum field theories are affected by UV divergences that invalidate 
simple algebraic proofs. 

The origin of UV divergences in field theory is double. First a field contains 
an infinite number of degrees of freedom. The corresponding divergences are 
directly related to renormalization group and reflect the property that even in 
renormalizable quantum fleld theories degrees of freedom remain coupled on all 
scales. 

However another of type of divergences can appear, which is related to the 
order between quantum operators and the transition between classical and quan- 
tum hamiltonians. Such divergences are already present in perturbation theory 
in ordinary quantum mechanics, for instance in the quantization of the geodesic 
motion on a manifold (like a sphere) . Even in the case of forces linear in the ve- 
locities (like a coupling to a magnetic fleld) flnite ambiguities are found. In local 
quantum fleld theory the problem is even more severe because for a scalar fleld 
for example the commutation between fleld operator (p and conjugate momentum 
TT takes the form 

[4>ix),7T{y)]=ihd''-\x-y). 

In a local theory all operators are taken at the same point and thus the commu- 
tator is divergent except in quantum mechanics [d = 1 with our conventions). 

Divergences of this nature thus are present as soon as derivative couplings are 
involved, or when fermions are present. They reflect the property that the knowl- 
edge of the classical theory is not sufficient, in general, to completely determine 
the quantized theory. 

Regularization. Regularization is a useful intermediate step in the renormal- 
ization program that consists in modifying the initial theory at short distance, 
large momentum or otherwise to render perturbation theory finite. Note that 
from the point of view of Particle Physics all these modifications alter in some es- 
sential way the physical properties of the theory, and thus can only be considered 
as intermediate steps in the removal of divergences. 

When a regularization can be found which preserves the symmetry of the initial 
classical action, a symmetric quantum field theory can be constructed. 

Momentum cut-off regularization schemes, based on modifying propagators at 
large momenta, are specifically designed to cut the infinite number of degrees 
of freedom. With some care these methods will preserve formal symmetries 
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of the un-renormalized theory that correspond to global (space- independent) 
linear group transformations. Problems may, however, arise when the symmetries 
correspond to non-linear or local transformations, like in the examples of non- 
linear o models or gauge theories, due to the unavoidable presence of derivative 
couplings. It is easy to verify that in this case regular izations that only cut 
momenta do not in general provide a complete regularization. 

The addition of regulator fields has in general the same effect as modifying 
propagators, but offers a few new possibilities, in particular when regulator fields 
have the wrong spin-statistics connection. Fermion loops in a gauge background 
can be regularized by such a method. 

Other methods have to be explored. In many examples dimensional regu- 
larization solves the problem because then the commutator between field and 
conjugated momentum taken at the same point vanishes. However in the case 
of chiral fermions dimensional regularization fails because chiral properties are 
specific to even dimensions. 

Of particular interest is the method of lattice regularization, because it can 
be used, beyond perturbation theory, either to discuss the existence of a quan- 
tum field theory, or to determine physical properties of field theories by non- 
perturbative numerical techniques. One verifies that such a regularization indeed 
specifies an order between quantum operators. It therefore solves the ordering 
problem in non-linear (j-models or non-abelian gauge theories. However, again 
it fails in presence of chiral fermions: the manifestation of this difficulty takes 
the form of a doubling of the fermion degrees of freedom. Until recently this had 
prevented a straightforward numerical study of chiral theories. 

Anomalies. That no conventional regularization scheme can be found in the 
case of gauge theories with chiral fermions is not surprising since we know exam- 
ples of theories with anomalies, i. e. theories in which a local symmetry present 
at the in the tree or classical approximation cannot be implemented in the full 
quantum theory. Note that this creates obstructions to the construction of chiral 
gauge theories because exact gauge symmetry, and thus the absence of anomaly, 
is essential for the physical consistency of a gauge theory. 

Note that we study in these lectures only local anomalies, which can be deter- 
mined by perturbative calculations; peculiar global non-perturbative anomalies 
have also been exhibited. 

Anomalies are local quantities because they are consequences of short distance 
singularities. They are responses to local (space-dependent) group transforma- 
tions but vanish for a class of space-independent transformations. This gives 
them a topological character that is further confirmed by their relations with the 
index of the Dirac operator in a gauge background. 

The recently discovered solutions of the Ginsparg-Wilson relation and the 
method of overlap fermions seem to provide an unconventional solution to the 
regularization problem in gauge theories with chiral fermions on the lattice. They 
evade the doubling problem of fermion because chiral transformations are no 
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longer strictly local on the lattice, and relate the problem of anomalies with 
the invariance of the fermion measure. The absence of anomalies can then be 
verified directly on the lattice, and this seems to confirm that the theories that 
had been discovered anomaly-free in perturbation theory are also anomaly-free 
in the non-perturbative lattice construction. Therefore the specific problem of 
lattice fermions was in essence technical rather than reflecting an inconsistency 
of chiral gauge theories beyond perturbation theory, as one may have feared. 

Finally since these new regularization schemes have a natural implementation 
in five dimensions under the form of domain wall fermions, it again opens the 
door to speculations about higher space dimensions. 

We first discuss the advantages and shortcomings, from the point of view of 
symmetries, of three regularization schemes, momentum cut-off, dimensional and 
lattice regularization. We show that they leave room for possible anomalies when 
both gauge fields and chiral fermions are present. 

We then recall the origin and the form of anomalies, beginning with the sim- 
plest example of the so-called abelian anomaly, i.e. the anomaly in the conserva- 
tion of the abelian axial current in gauge theories. We relate anomalies to the 
index of a covariant Dirac operator in the background of a gauge field. 

In the two-dimensional CP{N—1) models and in four-dimensional non-abelian 
gauge theories we exhibit gauge instantons. We show that they can be classi- 
fied in terms of a topological charge, space integral of the chiral anomaly. The 
existence of gauge field configurations that contribute to the anomaly has direct 
physical implications, like possible strong CP violation and the solution to the 
U{1) problem. 

We examine the form of the anomaly for a general axial current, and infer 
conditions for gauge theories that couple differently to fermion chiral components 
to be anomaly-free. A few physical applications are also briefly mentioned. 

Finally the formalism of overlap fermions on the lattice and the role of the 
Ginsparg-Wilson relation are explained. The alternative construction of domain 
wall fermions is explained, starting from the basic mechanism of zero-modes in 
super symmetric quantum mechanics. 

Conventions. Throughout these notes we work in euclidean space (with imag- 
inary or euclidean time), and this also implies a formalism of euclidean fermions. 

2 Momentum cut-off regularization 

We first discuss methods that work in the continuum (compared to lattice meth- 
ods) and at fixed dimension (unlike dimensional regularization). The idea then 
is to modify field propagators beyond a large momentum cut-off to render all 
Feynman diagrams convergent. However the regularization has to satisfy one 
important condition: the inverse of the regularized propagator must remain a 
smooth function of the momentum p. Indeed singularities in momentum vari- 
ables generate, after Fourier transformation, contributions to the large distance 
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behaviour of the propagator, and regularization should modify the theory only 
at short distance. 

2. 1 Matter fields: propagator modification 

Scalar fields. A simple modification of the propagator improves the convergence 
of Feynman diagrams at large momentum. For example in the case of the action 
of the self-coupled scalar field, 

S{ct>) = j d^x [i0(x)(-V2 +m2)0(x) + , (2.1) 
the propagator in Fourier space l/iw? -\- p^) can be replaced by 

with 

n 

Ab'(p) = (/ + m") n(l +pVMf). (2.2) 

4=1 

The masses are proportional to the momentum cut-off A, 

Mi = ctiA , cti > . 

If the degree n is chosen large enough all diagrams become convergent. In the for- 
mal large cut-off limit, at parameters a fixed, the initial propagator is recovered. 
This is the spirit of momentum or Pauli-Villars's regularization. 

Note that such a propagator cannot be derived from a hermitian hamiltonian. 
Hermiticity of the hamiltonian implies that if the propagator is, as above, a 
rational function, it must be a sum of poles with positive residues and thus 
cannot decrease faster than 1/p^. 

While this modification can be implemented also in Minkowski space because 
the regularized propagators decreases in all complex p'^ directions (except real 
negative), in euclidean time more general modifications are possible. Schwinger's 
proper time representation suggests: 

/■oo 

Ab{p)= dtp{tK^)e-*^p"+^"\ (2.3) 

^0 

in which the function p{t) is positive (to ensure that Ab(p) does not vanish and 
thus is invertible) and satisfies the condition 

|1 - p{t)\ < Ce-"^* ((7 > 0) for t^+oo. 
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By choosing a function p(t) that decreases fast enough for t ^ the behaviour 
of the propagator can be arbitrarily improved. If p{t) = 0{t^) the behaviour 
( |2.2| ) is recovered. Another example is: 



pit) = e{t-i), (2.4) 

9{t) being the step function, which leads to exponential decrease: 

e-(p2+m2)/A^ 

Ab(p) = 2^ 2 • (2-5) 

As the example shows, it is thus possible to find in this more general class prop- 
agators without unphysical singularities, but they do not follow from a hamilto- 
nian formalism because continuation to real time becomes impossible. 

Spin 1/2 fermions. For spin 1/2 fermions similar methods are applicable. To 
the free Dirac action 



5po = J d''xtp{x){0 + m)tp{x) , (2.6) 

corresponds in Fourier representation the propagator l/{m + i'^). We can replace 
it by the regularized propagator Af(p), 

n 

A-\p) = {m + z^) Hil + p^Mf). (2.7) 

i=l 

Note that we use the standard notation ^ = p^.'^^, with euclidean fermion 
conventions, analytic continuation to imaginary or euclidean time of the usual 
Minkowski fermions, and hermitian matrices 7^. 

Remarks. Momentum cut-off regularizations have several advantages: one 
can work at fixed dimension and in the continuum. However, two potential 
weaknesses have to be stressed: 

(i) The generating functional of correlation function Z{J) obtained by adding 
to the action f^Tl]) a source term for fields 



S{(t)) ^ S{(j)) - j d^x J{x)(t){x), 

can be written, 

Z(J) =deti/2(AB)exp[-Vi(V5J)]exp j A'^x d'^y J{x)A^{x - y)J{y) 



(2.8) 
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where the determinant is generated by the gaussian integration, and 



None of the momentum cut-off regularizations described so far can deal with the 
determinant. As long as the determinant is a divergent constant that cancels 
in normalized correlation functions this is not a problem, but in the case of a 
determinant in the background of an external field (which generates a set of 
one-loop diagrams) this may become a serious issue. 

(ii) This problem is related to another one: Some models have even in simple 
quantum mechanics divergences or ambiguities due to problem of order between 
quantum operators in products of position and momentum variables. A class 
of Feynman diagrams then cannot be regularized by this method. Quantum 
field theories where this problem occurs include models with non-linear or gauge 
symmetries. 

Global linear symmetries. To implement symmetries of the classical action in 
the quantum theory we need a regularization scheme that preserves the sym- 
metry. This requires some care but can always be achieved for linear global 
symmetries, i.e. symmetries that correspond to transformations of the fields of 
the form 

^fl(x) = R(t){x) , 

where is a constant matrix. The main reason is that in the quantum hamil- 
tonian field operators and conjugate momenta are not mixed by the transforma- 
tion, and therefore the order of operators is to some extent irrelevant. To take 
an example directly relevant here, a theory with massless fermions may, in four 
dimensions, have a chiral symmetry 



The substitution ( |2.7| ) (for m = 0) preserves chiral symmetry. Note the impor- 
tance here of being able to work in fixed dimension four because chiral symmetry 
is defined only in even dimensions. In particular the invariance of the integration 
measure [d1l'{x)d^l'{x)] relies on the property that tr75 = 0. 

2.2 Regulator fields 

Regularization in the form ( p.2| ) or ( |2.7| ) has another equivalent formulation 
based on the introduction of regulator fields. 



The scalar case. In the case of scalar fields to regularize the action ( |2.1|) for 
the scalar field 4) one introduces additional dynamical fields (f)r, r = 1, . . . , rmax, 
and considers the modified action S^egS't'^ 4'k)- 

5reg. {(P, (Pr) = j ^''x (- + m^) + ^ ^<P, (-V^ + M^) 

+VlicP+j:r<Pr)]. (2.9) 
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With the action (|2.9| ) any internal (p propagator is replaced by the sum of the (p 
propagator and all the (pr propagators Zr/{p'^ + M^). For an appropriate choice 
of the constants Zr, after integration over the regulator fields, the form ( |2.2| ) is 
recovered. Note that the condition of cancellation of the contribution at 
large momentum implies 

r 

Therefore not all Zr can be positive and some of the fields (pr necessarily are 
unphysical. 

Fermions. The fermion inverse propagator (|2.7|) can be written 



A-\p) = {m + z/) n (1 + i^/Mr){l - i^/Mr). 



r=l 



This indicates that again the same form can be obtained by a set of regulator 
fields {^r±i^r±}- One replaces the kinetic part of the action by 



(i'^xip{x){0 m)tp{x) ^ J d'^x 'ip{x){0 + 771)1/} (x) 

+ J2 — [ <i'^xi^reix)i0 + eMr)tpreix).i2.10) 

In the same way in the interaction the fields tp and i/j are replaced by the sums 

tp tp + ^'tpre, tp tp + ^tpre ■ 



Again for a proper choice of the constants Zr, after integration over the regulator 
fields the form ( |2.7| ) is recovered. Note in particular that for m = one finds 

Zj — |- = Zj- This indicates how chiral symmetry is preserved by the regulariza- 

tion, although the regulators are massive: by fermion doubling. The fermions 
ip+ and tp- are chiral partners. For a pair ip = (-0+, ?/'_), tp = {ipj^,ip-) the 
action can be written 



j d'^x ipix) (g) 1 + Ml (g) as) ipix) 



where the first matrix 1 and the Pauli matrix act in ± space. The spinors 
then transform like 



because cxi anticommutes with as. 
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2.3 Abelian gauge theory 

The problem of matter in presence of a gauge field can be decomposed into 
two steps, first matter in an external gauge field, and then the integration over 
the gauge field. For gauge fields we choose a covariant gauge, in such a way that 
power counting is the same as for scalar fields. 

Charged fermions in a gauge background. The new problem that arises in 
presence of a gauge field is that only covariant derivatives are allowed, because 
gauge invariance is essential for the physical consistency of the theory. The 
regularized action in a gauge background now reads 

S{iP,iP,A)= I d''x^{x){m+V))\[{l-V>-' /Ml)i^{x), (2.11) 

where is the covariant derivative 

Dm = + ^e^P • 

Note that up to this point the regularization, unlike dimensional or lattice reg- 
ularizations, preserves a possible chiral symmetry for m = 0. 

The higher order covariant derivatives however generate new, more singular, 
gauge interactions and it is no longer clear whether the theory can be rendered 
finite. 

Fermion correlation functions in the gauge background are generated by: 



Z{r],r];A)— / [dip{x)d'ip{x)] exp 



S{ijj,ijj, A) + I d xr]{x)ijj{x) + ^jJ{x)r]{x) 

(2.12) 

where rj, r] are Grassmann sources. Integrating over fermions explicitly we obtain 



Z{r],r];A) = Zo{A)exp 



- I d'^xd%r]{y)AF{A;y,x)ri{x) 



(2.13) 



Zo(A) ^TVdet 



(m+]Z))n(l-0V^') 



where A/" is a gauge field-independent normalization and A-f{A; y, x) the fermion 
propagator in an external gauge field. 

Diagrams constructed from Af(^; y, x) belong to loops with gauge field prop- 
agators, and therefore can be rendered finite if the gauge field propagator can 
be improved, a condition that we check below. The other problem involves the 
determinant that generates closed fermion loops in a gauge background. Using 
In det = tr In we find 



lnZo{A) = trln(m +P) + ^trln (l -pV^r) - = 0)> 

r 
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or using the anticommutation of 75 with p 

det(p + m) = det75(p + 771)75 = det(m -p), 
InZo(A) = i trln (m^ -p^) + ^trln (l -p^/^r) - = 0), 

r 

We see that the regularization has no effect from the point of view of power 
counting on the determinant, and therefore on one-loop diagrams of the form 
of fermion closed loops with external gauge fields, a problem that requires an 
additional regularization. 

The fermion determinant. The fermion determinant can finally be regularized 
by adding to the action a boson regulator field with fermion spin, and therefore 
a propagator similar to Ap but with different masses 

SBi^,<P;A) = I d''x^ix){Mi +]d)ll{l-V>yiM^f)Hx). (2.14) 

r=l 

The integration over the boson ghost fields (j), (j) adds to In Zq the quantity 
51nZo(A) = -itrln((Mo^)2-p2) _^trln(l-pV(M,B)2) -(A = 0). 

r=l 

Expanding in inverse powers of ]3 one adjusts the masses to cancel as many 
powers as possible. However, the unpaired initial fermion mass m is the source 
of a problem. The corresponding determinant can only be regularized with an 
unpaired boson . In the chiral limit m = we have two options: either we 
give a chiral charge to the boson field and the mass breaks chiral symmetry, 
or we leave it invariant in a chiral transformation. Then we find the determinant 
of the transformed operator 

For 9{x) constant e^^'^^ anticommutes with p and cancels. Otherwise a non- 
trivial contribution remains. The method thus suggests possible difficulties with 
space-dependent chiral transformations. 

Since actually the problem reduces to the study of a determinant in an ex- 
ternal background one can study it directly, as we will starting with section ^. 
One examines whether it is possible to define some regularized form in a way 
consistent with chiral symmetry. When this is possible one then inserts the one- 
loop renormalized diagrams in the general diagrams regularized by the preceding 
cut-off methods. 

Boson determinant in a gauge background. The boson determinant can be 
regularized by introducing a massive scalar charged fermion. It can also be 
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expressed in terms of the statistical operator using Schwinger's representation 
(trln = Indet) 

/"°° df 

In det H-ln det Hq = ti — [e'^^^ - e"*^] , 

Jo ^ 

where the operator H is analogous to a non-relativistic hamiltonian in a magnetic 
field, 

H = -D^D^+m^, Ho = ~V^+m^. 

The integral over time is regularized by cutting it for t small, integrating from 
t= 1/A2. 

The gauge propagator. For the free gauge action in a covariant gauge usual 
derivatives can be used because in an abelian theory the gauge field is neutral. 
The tensor F^i, is gauge invariant and the action for the scalar 9^^^ is arbitrary. 
Therefore the large momentum behaviour of the gauge field propagator can be 
arbitrarily improved. 

2.4 Nan- abelian gauge theories 

Compared with the abelian case, the new features of the non-abelian gauge 
action are the presence of gauge field self- inter actions and ghost terms. For 
future purpose wc define our notation. We introduce the covariant derivative, as 
acting on matter field, 

D^ = a^ + A^(a;), (2.15) 
where is an anti-hermitian matrix, and the curvature tensor F^j^ 

F^^ = [D^,D^] = d^A^ - d^A^ + [A^, A^]. (2.16) 

The pure gauge action then is 

S{A^) = -4^tr j d'^x trF^4x)F^4x), (2.17) 

In the covariant gauge 

the ghost field action takes the form: 

5ghost(A^, C,C) = - j d^x trCa^ {d^C + [A^, C]) . 

The ghost fields thus have a simple 5ab/p^ propagator and canonical dimension 
one in four dimensions. 
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The problem of regularization in non-abelian gauge theories has several fea- 
tures in common with the abelian CclSGj clS well as with the non-linear u-model. 
The regularized gauge action takes the form: 

5reg.(A^) = -/d'^^ trF^^P (D7 A2) F^, , (2.18) 

in which P is a polynomial of arbitrary degree. In the same way the gauge 
function 5^A^ is changed into: 

d^^^ ^ Q (^7 A2) d^A^ , (2.19) 

in which Q is a polynomial of same degree as P. As a consequence both the gauge 
field propagator and the ghost propagator can be arbitrarily improved. However, 
as in the abelian case, the covariant derivatives generate new interactions that are 
more singular. It is easy to verify that the power counting of one-loop diagrams 
is unchanged while higher order diagrams can be made convergent by taking the 
degrees of P and Q large enough: Regularization by higher derivatives takes care 
of all diagrams except, as in all geometric models, some one-loop diagrams (and 
thus subdiagrams) . 

As with charged matter the one-loop diagrams have to be examined separately. 
For fermion matter it is however still possible as, in the abelian case, to add a 
set of regulator fields, massive fermions and bosons with spin. In the chiral 
situation the problem of the compatibility between the gauge symmetry and the 
quantization is reduced to an explicit verification of the WT identities for the 
one-loop diagrams. Note that the preservation of gauge symmetry is necessary 
for the cancellation of unphysical states in physical amplitudes, and thus essential 
to the physical relevance of the quantum field theory. 

3 Other regularization schemes 

The other regularization schemes we now discuss have the common property 
that they modify in some essential way the structure of space-time, dimension 
regularization because it relies on defining Feynman diagrams for non-integer di- 
mensions, lattice regularization because continuum space is replaced by a discrete 
lattice. 

3. 1 Dimensional regularization 

Dimensional regularization involves continuation of Feynman diagrams in the 
parameter d {d is the space dimension) to arbitrary complex values, and therefore 
seems to have no meaning outside perturbation theory. However this regulariza- 
tion very often leads to the simplest perturbative calculations. 
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In addition it solves with the problem of commutation of quantum operators 
in local field theories. Indeed commutators for example in the case of a scalar 
field take the form 

[4>{x), n{y)] = ihS'^-^x -y) = ih{2T^f-'^ j d'^" V e^^^^-^^) , 

where ti{x) is the momentum conjugated to the field (l){x). As we have already 
stressed, in a local theory all fields are taken at same point, and therefore a 
commutation in the product (/)(a;)7r(a;) generates a divergent contribution (for 
d> 1) proportional to 

(5'^-i(0) = {2ttY-'^ j d'^-V- 

The rule of dimensional regularization imply the consistency of the change of 
variables p ^ \p and thus J d'^p = 0, in contrast to momentum regularization 
where it is proportional to a power of the cut-off. Therefore the order between 
operators becomes irrelevant because the commutator vanishes. Dimensional 
regularization thus is applicable to geometric models where these problems of 
quantization occur, like non-linear cr models or gauge theories. 

Its use, however, requires some care in massless theories. For instance in a 
massless theory in two dimensions integrals of the form J d^k/k'^ are met. They 
again vanish in dimensional regularization for the same reason. However here 
they correspond to an unwanted cancellation between UV and IR logarithmic 
divergences. 

More important here, it is not applicable when some essential property of 
the field theory is specific to the initial dimension. An example is provided by 
theories containing fermions in which parity symmetry is violated. 

Fermions. For fermions belonging to the fundamental representation of the 
spin group Spin((i) the strategy is the same. The spin problem can be reduced to 
the calculation of traces of 7 matrices. Therefore only one additional prescrip- 
tion for the trace of the unit matrix is needed. There is no natural continuation 
since odd and even dimensions behave differently. Since no algebraic manipu- 
lation depends on the explicit value of the trace, any smooth continuation in 
the neighbourhood of the relevant dimension is satisfactory. A convenient choice 
is to take the trace constant. In even dimension as long as only 7^ matrices 
are involved no other problem arises. However no dimensional continuation that 
preserves all properties of 7d+i can be found. This leads to serious difficulties if 
7d+i in the calculation of Feynman diagrams has to be expressed in terms of the 
product of all other 7 matrices. For example in four dimensions 75 is related to 
the other 7 matrices by 



4! 75 — . . ./i4 7mi ■ ■ ■ TAt4 : 



(3.1) 
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where e^^...^^ is the complete antisymmetric tensor with €1234 = 1. Problems 
therefore arise in the case of gauge theories with chiral fermions, because the 
special properties of 75 are involved as we recall below. This difficulty is the 
source of chiral anomalies. 

Since perturbation theory involves the calculation of traces, one possibility is 
to define 75 near four dimensions by 

75 = ^mi-M4 7mi---7m4, (3-2) 

where E^^p^ is a completely antisymmetric tensor, which reduces to — e^!ypcr/4! 
in four dimensions. It is easy to then verify that, with this definition, 75 anti- 
commutes with the other 7^ matrices only in four dimensions. If for example we 
evaluate the product ^vl^lv in d dimensions, we find: 

71^757^ = (f^ - 8)75 ■ 

Anticommuting properties of the 75 would have led to a factor —d instead. Mul- 
tiplied by a pole in l/((i— 4) consequence of UV divergences in one-loop diagrams, 
the additional contribution proportional to (i — 4 may lead to a finite difference 
with the formal result. 

The other option would be to keep the anticommuting property of 75, but the 
preceding example shows that this is contradictory with a form (|3.2|) . Actually 
one verifies that the only consistent prescription for generic dimensions then is 
that the traces of 75 with any product of 7^ matrices vanishes, and therefore is 
useless. 

Finally an alternative possibility consists in breaking 0{(£) symmetry and keep- 
ing the four 7 matrices of (i = 4. 

3.2 Lattice regularization 

We have explained that Pauli-Villars's regularization does not work for field 
theories in which the geometric properties generate interactions like models on 
homogeneous spaces (for example the non-linear cr-model) or gauge theories. In 
these theories some divergences are related to the problem of quantization and 
order of operators, which already appears in simple quantum mechanics. Other 
regularization methods then are needed. In many cases lattice regularization 
may be used. 

Lattice field theory. To each site x of the lattice are attached field variables 
corresponding to fields in the continuum. To the action S in the continuum 
corresponds a lattice action, the energy of lattice field configurations in the lan- 
guage of classical statistical physics. The functional integral becomes a sum over 
configurations and the regularized partition function is the partition function of 
a lattice model. 

All expressions in these notes will refer implicitly to a hypercubic lattice and 
we denote by a the lattice spacing. 
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The advantages of lattice regularization are: 

(i ) Lattice regularization indeed corresponds to a specific choice of quantiza- 
tion. 

(ii) It is the only established regularization that for gauge theories and other 
geometric models has a meaning outside perturbation theory. For instance the 
regularized functional integral can be calculated by numerical methods, like 
stochastic methods (Monte-Carlo type simulations) or strong coupling expan- 
sions. 

(iii) It preserves most global or local symmetries with the exception of the 
space 0{d) symmetry that is replaced by a hyper cubic symmetry (but this turns 
out not to be a major difficulty) and fermion chirality, which turns out to be a 
more serious problem, as we will show. 

The main disadvantage is that it leads to rather complicated perturbative 
calculations. 



3.3 Boson field theories 

Scalar fields. To the action (J2.1D for a scalar field (p in the continuum corresponds 



a lattice action, which is obtained in the following way: The euclidean lagrangian 
density becomes a function of lattice variables (p^x), where x now is a lattice site. 
Locality can be implemented by considering lattice lagrangians that depend only 
on a site and its neighbours (though this is a too strong requirement; lattice 
interactions decreasing exponentially with distance are also local). Derivatives 
dfj,(f) of the continuum are replaced by finite differences, for example: 

d^(j) ^ V;:'-(/) = [(t>{x + an^) - (/)(x)] /a , (3.3) 

where a is the lattice spacing and the unit vector in the fj, direction. The 
lattice action then is the sum over lattice sites. 

With the choice ( p.3| ) the propagator Aa{p) for the Fourier components of a 
massive scalar field is given by 

2 

A-\p) = + - 5^(1 - cos [ap^)). (3.4) 

It is a periodic function of the components of the momentum vector with 
period 2n/a. In the small lattice spacing limit the continuum propagator is 
recovered: 

ip)=m'+p'-^Yl ^'Pt + O {pD ■ (3.5) 

In particular hypercubic symmetry implies 0{d) symmetry at order p^. 

Gauge theories. Lattice regularization defines unambiguously a quantum the- 
ory. Therefore, once one has realized that gauge fields should be replaced by link 
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variables corresponding to parallel transport along links of the lattice, one can 
regularize a gauge theory. 

The link variables \Jxy are group elements associated with the links joining 
the sites x and y on the lattice. The regularized form of / dx is the product 
of link variables along a closed curve on the lattice, the simplest being a square 
on a hypercubic lattice, leading to the well-known plaquette action, each square 
forming a plaquette. The typical gauge invariant lattice action corresponding to 
the continuum action of a gauge field coupled to scalar bosons then has the form: 

plaquettes links sites 

(3.6) 

where x, y,... denotes lattice sites, and (3 and k are coupling constants. The ac- 
tion ( |3.6| ) is invariant under independent group transformations on each lattice 
site, lattice equivalents of the gauge transformations of the continuum theory. 
The measure of integration over the gauge variables is the group invariant mea- 
sure on each site. Note that on the lattice and in a finite volume the gauge 
invariant action leads to a well-defined partition function because the gauge 
group is compact. However in the continuum or infinite volume limits the com- 
pact character of the group is lost. Even on the lattice regularized perturbation 
theory is defined only after gauge fixing. 

We finally note that on the lattice the difficulties with the regularization do 
not come from the gauge field directly, but involve the gauge field only through 
the integration over chiral fermions. 

3.4 Fermion and the doubling problem 

We now review a few problems specific to relativistic fermions on the lattice. 
We consider the free action for a Dirac fermion: 

To regularize this action by a lattice and preserve chiral properties in the massless 
limit one can replace d^'ilj{x) by 

V|f*-'0(a;) = [i'{x + an^) - ip {x - an^)] /2a. 

Then the inverse of the fermion propagator A for the Fourier components i^{p) 
of the field is: 

A-\p) = m + ^J2^,'^ (3.7) 

^^ 

a periodic function of the components p^ of the momentum vector. A problem 
then arises: the equations relevant to the small lattice spacing limit. 



sin(ap^) = 
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have each two solutions = and = 7r/a within one period, i.e. within 
the BriUouin zone 2ti /a. Therefore the propagator ( p.7| ) propagates 2*^ fermions. 
To remove this degeneracy it is possible to add to the regularized action an 
additional scalar term 5S involving second derivatives: 

5S{i}, ip) = \M ^ [2t[){x)ip{x) -t[){x + an^) ip{x) - t[){x)ip {x + an^)] . (3.8) 

After Fourier transformation the modified Dirac operator Dw reads 

Dw{p) = m + (1 - cosap^) + - ^7^* sinap^ . (3.9) 

The fermion propagator becomes: 

Hp) = d\,{p) (dw{p)d\^{p)~^ 



with: 



Dw{p)DUp) 



m 



+ (1 - cos ap^) 



Therefore the degeneracy between the different states has been lifted. For each 
component that takes the value 7r/a the mass is increased by M. If M is of 
order 1/a the spurious states are eliminated in the continuum limit. This is the 
recipe of Wilson's fermions. 

However a problem arises if one wants to construct a theory with massless 
fermions and chiral symmetry. Chiral symmetry implies for the Dirac operator 
D{p) anticommutation with 75 

{i^(p),75} = 0, (3.10) 

and therefore both the mass term and the term (|3.8[) are excluded. It remains 
possible to add various counter-terms and try to adjust them to recover chiral 
symmetry in the continuum limit. But then there is no a priori guarantee that 
this is indeed possible and moreover calculations are plagued by fine tuning 
problems and cancellations of unnecessary UV divergences. 

One could also think about modifying the fermion propagator by adding terms 
connecting fermions separated by more than one lattice spacing. But it has been 
proven that this does not solve the doubling problem. In fact this doubling of 
the number of fermion degrees of freedom is directly related to the problem of 
anomalies. 

Since the most naive form of the propagator yields 2'^ fermion states, one tries 
in practical calculations to reduce this number to a smaller multiple of two, using 
for instance the idea of staggered fermions introduced by Kogut and Susskind. 

However the general picture has recently changed with the discovery of the 
properties of overlap fermions and solutions of the Ginsparg-Wilson relation, a 
topic we postpone, and on which we will come back in section |^. 
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4 The abelian anomaly 

We have pointed out that none of the standard regularization methods can deal in 
a straightforward way with one-loop diagrams in the case of gauge fields coupled 
to chiral fermions. We now show that indeed chiral symmetric gauge theories, 
involving gauge fields coupled to massless fermions, can be found where the axial 
current is not conserved. The divergence of the axial current in a chiral theory, 
when it does not vanish, is called an anomaly. Anomalies in particular lead to 
obstructions to the construction of gauge theories when the gauge field couples 
differently to the two fermion chiral components. 

Several examples are physically important like the theory of weak electro- 
magnetic interactions, the electromagnetic decay of the ttq meson, or the U{1) 
problem. We first discuss the abelian axial current, in four dimensions (the gen- 
eralization to all even dimensions then is straightforward), and then the general 
non-abelian situation. 

4.1 Abelian axial current and abelian vector gauge field 

The only possible source of anomalies are one-loop fermion diagrams in gauge 
theories when chiral properties are involved. This reduces the problem to the dis- 
cussion of fermions in background gauge fields, or equivalently to the properties 
of the determinant of the gauge covariant Dirac operator. 

We thus consider the QED-like fermion action S{ip, ip; A) for massless Dirac 
fermions ■0, in the background of an abelian gauge field A^ 

S{^,i};A) = - j d^x^{x)V)il}{x), Jd = + ie4{x), (4.1) 

and the corresponding functional integral 

Z{A^) = j [dil^di^] exp [-5(7/^, V^; A)] = detp . (4.2) 

We can find regularizations that preserve gauge invariance, 

^(x) = e^^(^)V'(a;), i^{x) = e-'^^''^'{x), A^{x) = --d^K{x) + A'^{x), 

(4.3) 

and since the fermions are massless, chiral symmetry. We would therefore naively 
expect the corresponding axial current to be conserved. However the proof of 
current conservation involves space-dependent chiral transformations, and there- 
fore steps that cannot be regularized without breaking the local symmetry. 
Under a space-dependent chiral transformation 



il:e{x) = e^^(^)T5 ij{x), Mx) = V'(^) e'^^^)^^ 



(4.4) 
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the action becomes 



where J?^{x), the coefficient of 9^6*, is the axial current, 



Jj(x) = itp{x)-f5-f^tp{x). 



(4.5) 



After the transformation (4^) Z^A^) becomes 



Z{A^,e) = det [e*T5^(^)pe^^5e(x) 



(4.6) 



Note that \n[Z{Afj,, 9)] is the generating functional of connected d^J^ correlation 
functions in an external field A^. 

Since e^'^^^ has a determinant that is unity, one would naively conclude that 
Z{A^,6) = Z{A^) and therefore that the current J^{x) is conserved. This is a 
conclusion we now check by an explicit calculation of the expectation value of 
dfj_J^{x) in the case of the action ([4.1|) . 

Remarks. 

(i) For any regularization that is consistent with the hermiticity of 75 



|Z(A^,^)r = det 



det 



e-*75e(x)pte-^T5^(^)l = det (pp 



and thus \Z{A^, 6)\ is independent of 6. Therefore an anomaly can appear only 
in the imaginary part of In 2. 

(ii) We have shown that one can find a regularization with regulator fields such 
that gauge invariance is maintained, and the determinant is independent of 6 for 
6{x) constant. 

(iii) If the regularization is gauge invariant Z{A^^6) is also gauge invariant. 
Therefore a possible anomaly will also be gauge invariant. 

(iv) In 2(^^,6*) is connected and IPI. Short distance singularities thus take 
the form of local polynomials in the fields and sources. Since a possible anomaly 
is a short distance effect it must also take the form of a local polynomial of A^ 
and d^9 constrained by parity and power counting, and 9^6* having dimension 
one and no mass being available, 

\nZ{A^,e) -\nZ{A^,Q) =i j d^x£(A, 9^;x), 

where C is the sum of monomials of dimension four. At order 6 only one is 
available: 

C{A,de]x) oc e^e^^p^d^e{x)A^{x)dpA^{x), 
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where e^i^p^ is the complete antisymmetric tensor with €1234 = 1. A simple 
integration by parts and anti-symmetrization shows that 

j A^xC{A, d9;x) (X e^e^^p^ j d'^x Fp^{x)Fp^{x), 

an expression that is gauge invariant. 

The coefficient of 9{x) is the expectation value in an external gauge field of 
dpjf^{x), the divergence of the axial current. It is thus determined up to a 
multiplicative constant, 

(dxJU^)) oc e'^e^^p^d^A^{x)dpA^{x) cx e^e^upaF^^{x)Fpa{x) , 

where Fp_j, = dp,Ai, — d^Ap, is the electromagnetic tensor, and we denote by (•) 
expectation values with respect to the measure e"*^*^^'^'"^). 

Since the possible anomaly is independent up to a multiplicative factor of the 
regularization, it must indeed be a gauge invariant local function of Ap. 

To find the multiplicative factor, which is the only regularization dependent 
feature, it is sufficient to calculate the coefficient of the term quadratic in A in 
the expansion of (9a in powers of A. We define the three-point function 

in momentum representation: 



6Ap{pi) SA^{p2] 
6 6 



A=0 

tr [757a]Z)-'(A;)] 



(4.7) 



SAp{pi) 5A^{p2 
r(3) is the sum of the two Feynman diagrams of figure |l| 



A=0 



k, Xm t 





(b) 



Fig. 1 Anomalous diagrams. 



■ Pi,/" 



The contribution of diagram (a) is: 



(a) 



(27r) 



tr 



and the contribution of diagram (b) is obtained by exchanging pi^ji, ^ P2t1u- 
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Power counting tells us that the function F*-^-* may have a linear divergence 
that, due to the presence of the 75 factor, must be proportional to e\^jjp, sym- 
metric in the exchange pi^'^v ^ P2, 7i/, and thus proportional to 

ex^,up{pi-P2)p ■ (4.9) 

On the other hand by commuting 75 we notice that T^^^ is formally a symmetric 
function of the three sets of external arguments. A divergence breaks the sym- 
metry between external arguments. Therefore a symmetric regularization of the 
kind we will adopt in the first calculation leads to a finite result. The result is 
not ambiguous because a possible ambiguity again is proportional to ([4.91). 

Similarly if the regularization is consistent with gauge invariance the vector 
current is conserved 

Pi^vf^Ak:PuP2) = Q. (4.10) 
Applied to the divergent part the equation implies 

-Pl^lP2p^x^llyp = , 

which cannot be satisfied. Therefore the sum of the two diagrams is finite. Finite 
ambiguities must also have the form ( |4.9| ) and thus are also forbidden by gauge 
invariance. All regularizations consistent with gauge invariance must give the 
same answer. 

Therefore there are two possibilities: 

(3) 

(i) The divergence k\r\^^{k;pi,p2) in a regularization respecting the sym- 
metry between the three arguments vanishes. Then both F*^'^-' is gauge invariant 
and the axial current is conserved. 

(ii) The divergence of the symmetric regularization does not vanish. Then it 
is possible to add to F*^^) a term proportional to (|4.9| ) to restore gauge invariance 
but this term breaks the symmetry between external momenta: the axial current 
is not conserved, an anomaly is present. 

4.2 Explicit calculation 

Momentum regularization. The calculation can be done using one of the various 
gauge invariant regularizations, for example momentum cut-off regularization or 
dimensional regularization with 75 being defined as in dimension four and thus no 
longer anticommuting with other 7 matrices. Instead we choose a regularization 
that preserves the symmetry between the three external arguments and global 
chiral symmetry, but breaks gauge invariance, modifying the fermion propagator: 

where e is the regularization parameter (e O-1-), p{z) is a positive differentiable 
function such that p(0) = 1, and decreasing fast enough for z — > +00, at least 
like 1/z. 
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Then current conservation and gauge invariance are compatible only if the 
divergence kxl'xij^j^ik; pi, P2) vanishes. 

It is convenient to consider directly the contribution C^'^\k) of order to 
{kxJ^{k)), which sums the two diagrams: 

X tr [iM+Hr'M^-^2r'iur'] , (4.11) 

because the calculation then suggests how the method generalizes to arbitrary 
even dimensions. 

We first transform the expression, using the identity: 



' + i^)-^ ^ 1 - ^{^ + (4.12) 
Then 

c^^\k) = ci^\k) + ci^\k), 

with 

Cl'\k) = e' ld'p,d'p2A^ip,)A,ip,) j ^^p{e{q + kf)p{e{q-p2f)p{eq^) 
xtr[757M(^-/2)-Vr'], (4.13) 



and 



j dVdV^^(pi)A.(p2) j ^Ae{q + kf)p{e{q-p2f)p{sq^) 

xti:[lM+^)~^lM-^2)-^l-4-^] ■ (4-14) 

In C2^{k) we use the cyclic property of the trace and the commutation of ^i,4~^ 
and 75 to cancel the propagator and obtain 

Ci'\k) = -e' J d''p,d''p2A^{p,)A,{p2) I 0^p{e{q + kf)p{e{q-p2f)p{eq^) 
xtr [757.(g^+|(')-SM(g^-/2)-'] . (4.15) 
We then shift q ^ q + P2 and interchange {pi, p) and (p2) 

Ci'\k) = -e' J d''p,d%AMMP2) I ^p{e{q-p2f)p{eq')p{e{q + p,f) 
xtr[757M(^-y2)-Vr'] • (4.16) 
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(2) (2) 

We see that the two terms and C2 would cancel in the absence of regulators. 
This corresponds to the formal proof of current conservation. However without 
regularization the integrals diverge and these manipulations are not legitimate. 

Instead here we find a non- vanishing sum due to the difference in regulating 
factors: 

C(^)^k)=e'l d'p^d'p2A,ip^)A^{p2) J -^p{e{q - P2f)p{eq') 

xtr [757^(3^-/2)- Vr'] [p{eiq + kf)-p{e{q+Pif)] • (4.17) 
After evaluation of the trace C*-^-* becomes (using ( p.l|) ): 

C^^\k) = -4e'J d^p^d''p2A,ip^)A^ip2) J -^p{e{q - p2r) p{eq') 

X e,.,.-^^^ [p{e{q + kf) - p{e{q + p^)')] . (4.18) 

Contributions coming from finite values of q cancel in the e —>■ limit. Due 
the cut-off the relevant values of q are of order We therefore rescale q 

accordingly qe^^"^ q and find 

C^^\k) = -Ae^ j d'p,d'p2A,{p,)A.{P2) j ^P((9-P2v^)')p(?') 

^ P2pqa P{{q + k^)^) - p{iq + PiVe)^) 

q^{q- P2\/eY ^fe 

Taking the £ ^ limit we obtain a finite result: 

C^'^\k) = -Ae'^e^.^p^ J d^pid^p2 A^(pi)A^(p2) V(pi, P2), 

with 

Ipa{pi,P2) ~ j j^^^^y^P2pqap'^ {eq'^) p' {eq'^) [2eqxik - pi) x] . (4.19) 
The identity: 

J dSq.qpfiq') = \S^p j d^q^fiq'), 
transforms the integral into: 

-±^p-'{eq^)p'{eq^). (4.20) 
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The remaining integral can be calculated explicitly (we recall p(0) = 1) 

and yields a result independent of the function p. We finally obtain: 

(kxJUk)) = -y^Cm-p- / dVd4p2Pi^A,(pi)p2pA,b2), (4.21) 
and therefore from the definition ( [4 .71) 

2 

^Ar\j^(/c;pi,P2) = ^e^,.pcTPipP2a ■ (4.22) 

This non- vanishing result implies that any definition of the determinant detp 
breaks at least either current conservation or gauge invariance. Since gauge 
invariance is essential to the consistency of a gauge theory we choose to break cur- 

(3) 

rent conservation. Exchanging arguments, we obtain the value of Pi^i^ xiiui^'i Pi^P'i) 

2 

PltJ^'^Xliuik'^Pl^P2) = -^^\upakpP2a ■ (4.23) 

If instead we had used a gauge invariant regularization, the result for V'^^^ would 
have differed by a term ^F^^^ proportional to 



5Tf^,{k-p,,p2)=Ktx^.p{pi -P2)p. (4.24) 
The constant K then is determined by the condition of gauge invariance 



Pn 

which yields 



0, 



2 

Pl^^^^fp..{k■puP2) = -^exupakpP2a ^ K = eV(67r2). (4.25) 
This gives an additional contribution to the divergence of the current 

k\5^XpAk]Pl,P2) = ^t^,XpaPlpP2a ■ (4.26) 

Therefore in a QED-like gauge invariant field theory with massless fermions the 
axial current is not conserved: this is called the chiral anomaly. For any gauge 
invariant regularization one finds 

k\'^flAk;Pl,P2) = = e^,^paPlpP2a ■ (4.27) 
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The equation ( [4.27|) can be rewritten after Fourier transformation as an axial 



current non-conservation equation: 

dxJlix) = -i—e^^p„F^^{x)Fp„{x) . (4.28) 

Since global chiral symmetry is not broken, the integral over the whole space 
of the anomalous term must vanish. This condition is indeed verified since the 
anomaly can immediately be written as a total derivative: 

^ IJ.V p<j fiv F pu Ad p^{tp^jy pijAjjd pAfj^ . (4.29) 

The space integral of the anomalous term depends only on the behaviour of 
the gauge field at boundaries, and this property already indicates a connection 
between topology and anomalies. 



The equation (|4.28|) also implies: 



In det 



e^-^^'^'^pe^^^^^-)] =lndetp-z^ j d^x e{x)t^,p,F^^{x)Fp^{x). (4.30) 



Remark. One might be surprised that in the calculation the divergence of the 
axial current does not vanish, though the regularization of the fermion propagator 
seems to be consistent with chiral symmetry. The reason is simple: if we add for 
example higher derivative terms to the action, the form of the axial current is 
modified and the additional contributions cancel the term we have found. 

In the form we have presented it the calculation generalizes without difficulty to 
general even dimensions 2n. Note simply that the permutation (pi, ii) ^ (p2, v) 
in equation ( |4.16D is replaced by a cyclic permutation. If gauge invariance is 
maintained the anomaly in the divergence of the axial current J fix) in general 
is 

s ^'^ 

d\J\ (x) = —2i ^^^^^^1 e^i!yi.../x„!/„-^^ii;yi • • -F^^jy^ , (4-31) 

where e^-^y-^...^„,y„ is the completely antisymmetric tensor, and = j^'^^^'> jg 
the axial current. 

Boson regulator. We have seen that we could also regularize by adding mas- 
sive fermions and bosons with spin, the unpaired boson affecting transformation 
properties under space-dependent chiral transformations. Denoting by the bo- 
son field and by M its mass, we perform in the regularized functional integral a 
change of variables of the form of a space-dependent chiral transformation acting 
in the same way on the fermion and boson field. The variation 5S of the action 
at first order in 9 is 

5S= I d^x [d^e{x)Jl{x) + 2zM^(x)^(x)75(^(x)] , 
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with 

J^{x) = itjj{x)'y5'yf,tp{x) + Z(/)(x)757^^(x). 

Expanding in 9 and identifying the coefficient of 9{x) we therefore obtain the 
equation 

{d^Jl{x)) = 2iM {^{x)j5(j){x)) = -2zMtr75 (x|p"^ |x) . (4.32) 

The divergence of the axial current comes here from the boson contribution. We 
know that in the large M limit it becomes quadratic in A. Expanding the r.h.s. in 
powers of A, keeping the quadratic term we find after Fourier transformation: 

C^'\k) = -2iMe^ I d^p^d^p2A^ip^)A,{p,) J ^ 

X tr [75(^+ ^ - ^M)-^-f^{^-^2 - tM)-^-f,{^ - iM)-^] . (4.33) 

The apparent divergence of this contribution is regularized by formally vanishing 
diagrams that we do not write, but which justify the following formal manipu- 
lations. 

In the trace the formal divergences cancel and one obtains 

C^^\k) ~ SM^eh^^p^ [d^Pid''p2PlpP2aAMMP2) 
1 [ d^q 



(27r)4 J (g2 + M2)3- 

The limit M — > oo corresponds to remove the regulator. The limit is finite 
because after rescaling of q the mass can be eliminated. One finds: 

C^'^\k) ~ — ^e^i^pcT / d^i dV {P2)pipP2aAp{pi)A^{p2) , 

in agreement with equation ( f4.27|) . 

Point- splitting regularization. Another calculation, based on regularization 
by point splitting, gives further insight into the mechanism that generates the 
anomaly. We thus consider the non-local operator 



J^{x, a) = iijj{x — 0/2)757^?/' (x -|- a/2) exp 



px+a/2 

ie / AA(s)dsA 

Jx-a/2 



'x-a/2 



(4.34) 



in the limit |a| — > 0. To avoid a breaking of rotation symmetry by the regulariza- 
tion, before taking the limit |a| ^ we will average over all orientations of the 
vector a. The multiplicative gauge factor (parallel transporter) ensures gauge 
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invariance of the regularized operator (transformations ( [4.3|) ). The divergence of 
the operator for |a| — > then becomes 



d^J^{x, a) ~ -eaxtpix - a/2)'y5'y^F^xix)tpix + a/2) 

f.x+a/2 



X exp 



ie / AA(s)dsA 

Jx-a/2 



(4.35) 



where the t/j, field equations have been used. We now expand the expectation 
value of the equation in powers of A. The first term vanishes. The second term 
is quadratic in A and yields 

{dlJl{x,a)) ^ ie''axF^,x{x) j d'^y A,{y+x)li^^A^{y-a/2)^,A^{-y-a/2)^^ , 
where Apd;) is the fermion propagator: 

We now take the trace. The propagator is singular for \y \ = 0{\a\) and therefore 
we can expand Ai,{x + y) in powers of y. The first term vanishes for symmetry 
reasons {y i— > — y), and we obtain 



{d^J^{x,a)) -^e^^raaxF^xix)dpA^{x) / d^y- 



y + a/2\^\y - a/2\^ 



The integral over y gives a linear combination of Sp^ and ttptta but the second 
term gives a vanishing contribution due to e symbol. It follows 

{dp',ix,a)) ~ ^,e,^r,axarF,x{x)dpAM J ^'^^TWlf^^W ' 
After integration we then find 

{d^Jlix^a)) ~ -^ef,^rp^^F^x{x)Fp^{x). 

Averaging over the a directions we see that the divergence is finite for \a\ 0, 
and thus finally 

■ 2 



l„l^o\ ' V 167r 

and we recover the result ( f4.28|) 
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On the lattice an averaging over is produced by summing over all lattice 
directions. Because the only expression quadratic in that has the symmetry 
of the lattice is the same result is found: the anomaly is lattice-independent. 

A direct physical application. In a phenomenological model of Strong Interac- 
tion physics, where the a SU{2) x SU{2) chiral symmetry is softly broken by the 
pion mass, in the absence of anomalies the divergence of the neutral axial cur- 
rent is proportional to the ttq field (corresponding to the neutral pion). A short 
calculation then shows the decay of ttq into two photons would vanish a zero mo- 
mentum. The axial anomaly ( [4.28| ) gives instead a non-vanishing contribution 



to the decay, in good agreement with experimental data. 

Chiral gauge theory. A gauge theory is consistent only if the gauge field is 
coupled to a conserved current. An anomaly that affects the current destroys 
gauge invariance at the quantum level. Therefore the theory with axial gauge 
symmetry, where the action in the fermion sector reads 

S{i;,i;;B) = - j d^^x ^{x){^ + ig-fMi.^). (4-37) 

is inconsistent. Indeed current conservation applies to the BBB vertex at one- 
loop order. Because now the three point vertex is symmetric the divergence is 
given by the expression ( [4.21] ), and thus does not vanish. 



More generally the anomaly prevents the construction of a theory that would 
have both an abelian gauge vector and axial symmetry, where the action in the 
fermion sector would read 

S{i),ip;A,B) = - j d'^xi}{x){0 + ie^ + i'-i^g$)ij{x). (4.38) 

A way to solve both problems is to cancel the anomaly by introducing another 
fermion of opposite chiral coupling. With more fermions other combinations of 
couplings are possible. Note, however that in the purely axial gauge theory it is 
easy to verify that a theory with two fermions of opposite chiral charges can be 
rewritten as a vector theory by combining differently the chiral components of 
both fermions. 

4.3 Two dimensions 

As an exercise we verify by explicit calculation the general expression ( [4. 311 ) in 
the special example of dimension two 

d^^Jl = -^^e^.F^, ■ (4.39) 



The general form of the l.h.s. is again dictated by locality and power counting: the 
anomaly must have canonical dimension two. The explicit calculation requires 
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some care since massless fields may lead to IR divergences in two dimensions. One 
thus gives a mass m to fermions, which breaks chiral symmetry explicitly, and 
takes the massless limit at the end of the calculation. The calculation involves 
only one diagram 



^=0 



(27r)2 ^ -^i^ + m'i^ + ii^ + m' 

Here the 7-matrices are simply the ordinary Pauli matrices. Then 

We use the method of the boson regulator, which leads to the two-dimensional 
analogue of equation ( |4.32| ). It here leads to the calculation of the difference 
between two diagrams (analogues of equation ( [4.33|) ) due to the explicit chiral 
symmetry breaking 



1 1 



C^{k) = 2m— — tr73 / d^g-^— -7^ , , - {m ^ M) 



In the trace again the divergent terms cancel 

The two contributions are now separately convergent. When m ^ the rm? 
factor dominates the logarithmic IR divergence and the contribution vanishes. 
In the second term in the limit M 00 one obtains 

C/i(^)L^o,M^oo ~ -4eM e^j^kj, ^^^^^ j ^^2 _^ M'^Y ^ -—(^i^fku , 
in agreement with equation ( f4.39| ). 

4.4 Non-abelian vector gauge fields and abelian axial current 

We still consider an abelian axial current but now in the framework of a non- 
abelian gauge theory. The fermion fields transform non-trivially under a gauge 
group G and is the corresponding gauge field. The action is: 

S{ip, tjj;A) = - j d^x 4^{x)^tP{x), (4.40) 
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with the convention ( p.l5|) and: 

P^0 + A- (4.41) 

In a gauge transformation of unitary matrix g(x) the gauge field and the 
Dirac operator become 

A^(a;) ^ gix)d^g-\x) + g{x)A^{x)g-\x) ^ p ^ g-\x)pg{x) . (4.42) 

The axial current J,^ 

J^(x) = zV'(x)757/,i/'(a^), 

is still gauge invariant. Therefore no new calculation is needed; the result is com- 
pletely determined by dimensional analysis, gauge invariance and the preceding 
abelian calculation that yields the term of order A^, 

dxJxix) = -z^e^^p^irY^^Yp^ , (4.43) 

in which F^,^ now is the corresponding curvature ( |2.16|) . Again this expression 
must be a total derivative. One indeed verifies: 

e^^p^ tr F^^Fp^ = 4 e^^pad^, tr( A^^p A^ + § A^ Ap A^) . (4.44) 



4-5 Anomaly and eigenvalues of the Dime operator 

We assume that the spectrum of p, the Dirac operator in a non-abelian gauge 
field (equation ( ^.41| )), is discrete (putting temporarily the fermions in a box if 
necessary) and call and ^nix) the corresponding eigenvalues and eigenvectors: 

= dn^n ■ (4.45) 

For a unitary or orthogonal group the massless Dirac operator is anti-hermitian; 
therefore the eigenvalues are imaginary and the eigenvectors orthogonal. In 
addition we choose them with unit norm. 

The eigenvalues are gauge invariant, because in a gauge transformation of 
unitary matrix g{x) the Dirac operator transforms like in ( [1.42|) , and thus simply 



Vn(a:) g{x)ipn{x). 

The anticommutation P75 + 75P = implies 

V>lb^n, = -dnl5V^n ■ (4.46) 

Therefore either dn is different from zero, and 75<^n is an eigenvector of p with 
eigenvalue —dn, or dn vanishes. The eigenspace corresponding to the eigenvalue 
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then is invariant under 75, which can be diagonahzed: the eigenvectors of can 
be chosen eigenvectors of definite chirahty, i.e. eigenvectors of 75 with eigenvalue 
±1, 

We call 7i_|_ and 7i_ the dimensions of the eigenspace of positive and negative 
chirality respectively. 

We now consider the determinant of the operator JZ) + m regularized by mode 
truncation (mode regularization) : 

det^vO + m) = JJ (rf^ + m), (4.47) 

n<Ar 

keeping the lowest eigenvalues of (in modulus), with — n+ — ri- even, in 
such a way that the corresponding subspace is 75 invariant. 

The regularization is gauge invariant because the eigenvalues of p are gauge 
invariant. 

Note that in the truncated space 

tr75 = n+-n_. (4.48) 

The trace of 75 equals n_(_ — n_, the index of the Dirac operator p. A non- 
vanishing index thus endangers axial current conservation. 

In a chiral transformation (|4.4| ) with constant 9 the determinant of (!0 + m) 
becomes 

detN(0 + m)^ detN {e'^'^' (0 + m) e*^'^^) . 

We now consider the various eigenspaces. 

If dn the matrix 75 is represented by the Pauli matrix ai in the sum of 
eigenspaces corresponding to the two eigenvalues ±dn and p + m by d^os + m. 
The determinant in the subspace then is 

det(e*^'"^(rf„a3 + m) e*^<"i) = dete^^^-^^ det((i^a3 + m) = - rf^, 

because ai is traceless. 

In the eigenspace of dimension uj^ of vanishing eigenvalues dn with eigenvectors 
with positive chirality, 75 is diagonal with eigenvalue 1 and thus 

m"+ ^ e2*^"+ . 
Similarly in the eigenspace of chirality —1 and dimension n_ 

We conclude 

detAr(e^^^^(P + m) e^^^^) ^ ^2^e{n+-u_) det;v(P + m). 
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The ratio of both determinants is independent of A^. Taking the Umit N ^ oo 
we find 



det 



(e^^5^(P + m) e^^s^) (0 + 



m] 



^2i6(n^ — n_ ) 



(4.49) 



Note that the l.h.s. of equation (|4.49[) is obviously 1 when 9 = nn, which imphes 
that the coefficient of 29 in the r.h.s. must indeed be an integer. 
The variation of In det (p + m) 



In det 



m] 



2i9 (n+ — n_) 



at first order in 9 is related to the variation of the action (|4.1J ) and thus to 
the expectation value of the integral of the divergence of the axial current, 
/ d^x(9^J^(x)) in four dimensions. In the limit m = it is thus related to 
the space integral of the chiral anomaly (|4.43| ) . 

We have thus found a local expression giving the index of the Dirac operator: 



327r 



2 '-fiiypcr 



d^x trF.^^F 



per 



n+ — n_ 



(4.50) 



Concerning this result several comments can be made: 

(i) At first order in 9 in the absence of regularization we have calculated 
(In det = tr In) 

In det [l + z^(75 + + m)750 + m)"^)] ~2z^tr75, 

where we have used the cyclic property of the trace. Since the trace of the 
matrix 75 vanishes we could expect naively a vanishing result. But trace here 
means trace in 7 space and in coordinate space, and 75 really stands here for 
755(x — y). The mode regularization yields a well-defined finite result for the ill 
defined product x 5^(0). 

(ii) The property that the integral ( ^.50| ) is quantized shows that the form 
of the anomaly is related to topological properties of the gauge field since the 
integral does not change when the gauge field is deformed continuously. The 
integral of the anomaly over the whole space thus depends only on the behaviour 
at large distances of the curvature tensor F„,^ and the anomaly must a total 



derivative as equation (|4.44 ) confirms. 

(iii) One might be surprised that det!0 is not invariant under global chiral 
transformations. However we have just established that when the integral of the 
anomaly does not vanish, det vanishes. This explains that to give a meaning to 



the r.h.s. of equation (|4.49|) we have been forced to introduce a mass to find a non- 
trivial result. The determinant of in the subspace orthogonal to eigenvectors 
with vanishing eigenvalue, even in presence of a mass, is chiral invariant by parity 
doubling, but for n+ ^ n_ not the determinant in the eigenspace of eigenvalue 
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zero because the trace of 75 does not vanish in the eigenspace (equation (|4.48|) ). 
In the hmit m —>■ the complete determinant vanishes but not the ratio of 
determinants for different values of 9 because the powers of m cancel. 

(iv) The discussion of the index of the Dirac operator is valid in any even 
dimension. Therefore the topological character and the quantization of the space 
integral of the anomaly are general. 

5 Instantons, anomalies and ^-vacua 

We now discuss the role of instantons in several examples where the classical 
potential has a periodic structure with an infinite set of degenerate minima. We 
exhibit their topological character, and in the presence of gauge fields relate 
them to anomalies and the index of the Dirac operator. Instantons imply that 
the eigenstates of the hamiltonian depend on an angle 9. In the quantum field 
theory the notion of ^-vacuum emerges. 

5.1 The periodic cosine potential 

As a first example of the role of instantons when topology is involved we consider 
a simple hamiltonian with a periodic potential: 

H = -^{d/dxf + ^sin'^x. (5.1) 

The potential has an infinite number of degenerate minima x = nTr, n G Z. 
Each minimum is an equivalent starting point for a perturbative calculation of 
the eigenvalues of H. Periodicity implies that the perturbative expansions are 
identical to all orders in g, a property that seems to imply that the quantum 
hamiltonian has an infinite number of degenerate eigenstates. In reality we know 
that the exact spectrum of the hamiltonian H is not degenerate, as a result of 
barrier penetration. Instead it is continuous and has, at least for g small enough, 
a band structure. 

The structure of the ground state. To characterize more precisely the struc- 
ture of the spectrum of the hamiltonian (|5.1|) we introduce the operator T that 
generates an elementary translation of one period tt 

T'il}{x) = ijjix + tt). 

Since T commutes with the hamiltonian, 

[T,if]=0, (5.2) 

both operators can be diagonalized simultaneously. Because the eigenfunctions of 
H must be bounded at infinity, the eigenvalues of T are pure phases. Each eigen- 
function of H thus is characterized by an angle 9 (pseudo-momentum) eigenvalue 
of T: 

T \9) ^ \9) . (5.3) 
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The corresponding eigenvalues En{6) are periodic functions of 6 and for g ^ Q 
are close to the eigenvalues of the harmonic oscillator 

EM =n + 1/2 + 0{g). 

To all orders in powers of g En{9) is independent of 9 and the spectrum of H is 
infinitely degenerate. Exponentially small contributions due to instantons lift the 
degeneracy and introduce a 9 dependence. To each value of n then corresponds 
a band when 9 varies in [0, 27r]. 

Path integral representation. The spectrum of H can be extracted from the 
calculation of the quantity Zi 



Z,i(3)=trT'e-^^ = l-f^ [d9 

^ n=0 



Indeed 



= J]e^^^Z,(/3) = J]e- 



■/3E„(9) 



(5.4) 



where Z{9, (3) is the partition function restricted to states with a fixed 9 angle. 

The path integral representation of Zi[l3) differs from the representation of 
the partition function Zq{I3) only by the boundary conditions. The operator T 
has the effect of translating the argument x in the matrix element {x' \ tre~^^ \x) 
before taking the trace: 



Ziil3) = / [dx(t)]exp[-5(a;)] 

J x{f3/2)=x{-l3/2)+e7T 

S{x) = — [x'^it) + sin^ (x(t))] dt . 

^9 J -13/2 



(5.5) 
(5.6) 



A careful study of the trace operation in the case of periodic potentials shows 
that x(— /9/2) varies over only one period (see Appendix |Xl]) . 

Therefore from (|5.4|) we derive the path integral representation of Z{9,(3): 



J X 



^ Jx{P/2)=x{-P/2)+liT 



[dx{i)\ exp [—S{x) + iL 



x(/3/2)=x(-/3/2) (mod tt) 



[dx(t)] exp 



9 f^/^ 
-S{x) +i- / dtx{t) 

^ J -13/2 



(5.7) 



Note that £ is a topological number since two trajectories with different values 
of £ cannot be related continuously. In the same way 



Q = - dtx(t), 

J-p/2 
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is a topological charge; it depends on the trajectory only through the boundary 
conditions. 

For P large and g the path integral is dominated by the constant solutions 
Xc{t) = mod n corresponding to the i = sector. A non-trivial 9 dependence 
can come only from instanton (non-constant finite action saddle points) contri- 
butions corresponding to quantum tunnelling. Note that quite generally 

j^[Mt)±Mmr>0 - 5> IcosW+oo)) -c„.(.(-.o))|A. (5.9) 

The action is finite only if x(±oo) G {0, n}. The non- vanishing value of the l.h.s. 
is 2. The minimum is reached for trajectories Xc solution of 

Xc = ±smxc ^ = 2tan-^e±(*-*o), (5.10) 

and the corresponding classical action then is: 

Six,) = 2/g. (5.11) 

The instanton solutions belong to the i = ±1 sector and connect two consec- 
utive minima of the potential. They yield the leading contribution to barrier 
penetration for g ^ 0. An explicit calculation yields 

Eo{g) = i?pert.(^7) - ^ cos^[l + 0{g)]. 



5.2 Instantons and anomaly: CP{N — 1) models 

We now consider field theories, the two-dimensional CP{N — 1) models, where 
again instantons and topology play a role and the semi- classical vacuum has a 
similar periodic structure. The new feature is the relation between the topolog- 
ical charge and the two-dimensional chiral anomaly. 

We here mainly describe the nature of the instanton solutions and refer the 
reader to the literature for a more detailed analysis. Note that the explicit cal- 
culation of instanton contributions in the small coupling limit in the CP{N — 1) 
models, as well as in the non-abelian gauge theories we discuss in section p.3| , 
remains to large extent an unsolved problem. Due to the scale invariance of 
the classical theory, instantons depend on a scale (or size) parameter. Instan- 
ton contributions then involves the running coupling constant at the instanton 
size. Both families of theories are UV asymptotically free. Therefore the run- 
ning coupling is small for small instantons and the semi-classical approximation 
is justified. However, in the absence of any IR cut-off, the running coupling 
becomes large for large instantons, and it is unclear whether a semi-classical 
approximation remains valid. 
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The CP{N — 1) manifolds. We consider a A^-component complex vector (f of 
unit length, 

(^■ip = l. (5.12) 

This condition characterizes a space isomorphic to the quotient space U{N) /U{N—^ 
1). In addition two vectors (p and (p' are considered equivalent if 

ip' = ip^ ip'^=e'^ip^. (5.13) 

This condition characterizes the symmetric space and complex Grassmannian 
manifold U{N) /U{1)/U{N - 1). It is isomorphic to the manifold CP{N - 1) 
(for — 1-dimensional Complex Projective), which is obtained from C"^ by the 
equivalence relation 

/ • f- / \ 

Zqi = Zq, II '^Zq^ 

where A belongs to the Riemann sphere (compactified complex plane). 

The CP{N — 1) models. A symmetric space admits a unique metric, up to a 
multiplicative factor, and this leads to a unique action with two derivatives. One 
form of the unique symmetric classical action is: 

S{<p,A^) = ^ j d^xD;;^-D^^, (5.14) 

in which is a coupling constant and the covariant derivative: 

B^ = d^ + iA^. (5.15) 

The field is a gauge field for the U{1) transformations 

^'(x) = e^^(^) ^{x) , A'^ix) = A^{x) - d^A{x). (5.16) 

The action is obviously U{N) symmetric and the gauge symmetry ensures the 
equivalence ( |5.13| ). 

Since the action contains no kinetic term for the gauge field is not a dy- 
namical field but only an auxiliary field that can be integrated out. The action 
is quadratic in A and the gaussian integration results in replacing in the action 
A^ by the solution of the A-field equation: 

A^ = iip-d^ip, (5.17) 

where the equation ( ^.12| ) has been used. After this substitution the field ■ d^(p 
acts as a composite gauge field. 

For what follows however we find more convenient to keep A^ as an indepen- 
dent field. 
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Instantons. To prove the existence of locally stable non-trivial minima of the 
action the following Bogomolnyi inequality can be used (note the analogy with 

j d^a; \D^^ T ie^^.'DM^ > , (5.18) 

(e^iy being the antisymmetric tensor, ei2 = 1). After expansion the inequality 
can be cast in the form 

S{ip)>27T\Q{ip)\/g, (5.19) 

with 

= -77-^^!^ / d^xD^if -D^if = ^ / d^xef^„D„Dfj_ip ■ If . (5.20) 
2n J in J 

Then 

ze^^D^D^ = ize^, [D,, D^] = \F^, , (5.21) 
where F^^ is the curvature 

Therefore, using ( |5.12| ), 

1 ^2. 



Q{^) = ^J d'xe^^F^,. (5.22) 

The integrand is proportional to the two-dimensional abelian chiral anomaly 
(|4.39| ) , and thus is a total divergence 

2^ -^fJ,!^ d^JL^^luA-ij . 



Substituting this form into equation (|5.22| ) and integrating in a large disk of 
radius R one obtains: 

Q{(p) = ^ lim (f dXf,A^{x). (5.23) 

Q((p) thus depends only on the behaviour of the classical solution for \x\ large and 
is a topological charge. Finiteness of the action demands that at large distances 
D^(/j vanishes and therefore 

Since 95 7^ this equation implies that F^^, vanishes and thus is a pure gauge 
(and (fi a gauge transform of a constant vector) 

A^ = a^A(x) Q{ip) = 7^ lim / dx^a^A(x) . (5.24) 

27r Ji^oo 
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The topological charge measures the variation of the angle A(x) on a large circle, 
which is a multiple of 27r because is regular. One is thus led to the consider- 
ation of the homotopy classes of mappings from U{1), i.e. Si to Si, which are 
characterized by an integer n, the winding number. This is equivalent to the 
statement that the homotopy group 7ri(S'i) is isomorphic to the additive group 
of integers Z. 
Then 

Q{ip) =n ^ S{ip) >2n\n\/g. (5.25) 

The equality S{ifi) = 2n\n\/g corresponds to a local minimum and implies that 
the classical solutions satisfy first order partial differential (self-duality) equa- 
tions: 

D^ip = ±ief^^D^(f . (5.26) 

For each sign there is really only one equation for instance = 1, z/ = 2. It is 
simple to verify that both equations imply the (^-field equations, and combined 
with the constraint ( [5.12| ) the A- field equation ( |5.17D . In the complex coordinates 
z — xi + ix2, z — xi — ix2 they can be written 

d:,(Pa{z,z) = -iAz{z,z)(po,{z,z), 
dz(Pa{z,z) = -iA^{z,z)ipa{z,z). 

Exchanging the two equations just amounts to exchange if and if. We therefore 
solve only the second equation 

iPc{z,z) = K{z,z)Pa{z), 

where k{z, z) is a particular solution of 

dzK{z,z) = —iAz{z,z)K,{z,z). 

Vectors solutions of the equations (|5.26| ) are proportional to holomorphic or 
anti-holomorphic (depending on the sign) vectors (this reflects the conformal 
invariance of the classical field theory). The function k{z,z), which gauge in- 
variance allows to choose real (this corresponds to the d^A^ = gauge), then is 
constrained by the condition ( p.l2| ) 

K^{Z,Z)P-P = 1. 

The asymptotic conditions constrain the functions Pa{z) to be polynomials. 
Common roots to all P^ would correspond to non-integrable singularities for ipa., 
and therefore are excluded by the condition of finiteness of the action. Finally if 
the polynomials have maximal degree n, asymptotically 

P^{z) ~ c«z" ^ ~ -^{z/-zr'\ 

Vc ■ c 



40 



When the phase of z varies by 27t, the phase of (/Pq varies by 2?7.7r, showing that 
the corresponding winding number is n. 

The structure of the semi- classical vacuum. In contrast to our analysis of 
periodic potentials in quantum mechanics, we have here discussed the existence 
of instantons without reference to the structure of the classical vacuum. To 
find an interpretation of instantons in gauge theories, it is useful to express the 
results in the temporal gauge. Then classical minima of the potential correspond 
to fields ip{xi), where xi is only the space variable, gauge transforms of a constant 
vector: 

=e^^(^i) V, v-v = l. 

Moreover if the vacuum state is invariant under space refiection (f{+oo) = 
(p{—oo) and thus 

A(+oo) - A(-oo) = 2z/7r z/ G Z . 

Again is a topological number that classifies degenerate classical minima, and 
the semi-classical vacuum thus has a periodic structure. This analysis is consis- 
tent with Gauss's law that only implies that states are invariant under infinites- 
imal gauge transformations, and therefore under gauge transformations of the 
class z/ = that are continuously connected to the identity. 

We now consider a large rectangle with extension R in the space direction and 
T in the euclidean time direction and by a smooth gauge transformation continue 
the instanton solution to the temporal gauge. Then the variation of the pure 
gauge comes entirely from the sides at fixed time. One finds for ^ oo, 

A(+oo, 0) - A(-oo, 0) - [A(+oo, T) - A(-oo, T)] = 2n7r . 

Therefore instantons interpolate between different classical minima. Like in the 
case of the cosine potential, one projects onto a proper quantum eigenstate, 
the "^-vacuum" corresponding to an angle 9 by adding, in analogy with the 
expression ( p.7|) , a topological term to the classical action 

9 f 

S{(p) S{ip) + i— / d^xe^jyF^^. 



Remark. Replacing in the topological charge Q the gauge field by the explicit 
expression ( p.l7| ) we find 

Qif) = 77- / e^^d^if ■ d^if = / d^x d(^a A dy?a , 
2n J Ztc J 

where the notation of exterior differential calculus has been used. We recognize 
the integral of a two-form, a symplectic form and 47rQ is the area of a 2-surface 
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embedded in CP{N — 1). A symplectic form is always dosed, here it is also 
exact, so that Q is the integral of a one-form (equation (|5.23| )) 



The 0(3) non-linear a-model. The CP(1) model is locally isomorphic to the 
0(3) non-linear u-model, with the identification 

f{x)=^^{x)al,pipfiix). (5.27) 

where cr* are the three Pauli matrices. 

Using for example an explicit representation of Pauli matrices, one verifies 
indeed 



fix)(P\x) = 1, d^(P\x)d^fix)=4B^ifB^ip. 

Therefore the field theory can be expressed in terms of the field and takes 
the form of the non-linear cr-model. The fields (p gauge invariant, and the 
whole physical picture is a picture of confinement of the charged scalar "quarks" 
(fia{x) and the propagation of the (f)^ neutral bound states. 

Instantons in the (p description take the form of (p configurations with uniform 
limit for \x\ oo. They thus define a mapping from the compactified plane 
topologically equivalent to 5*2 to the sphere 5*2 (the (p^ configurations). Since 
712 (^'2) = ^ the (p and (p pictures are consistent. 

In the example of CP{1) a solution of winding number 1 is 

1 z 



a/1 + ZZ y/T+ZZ 

Translating the CP(1) minimal solution into the 0(3) a-model language one 
finds 

z + z 1 z — Z 1 — ZZ 

(Pl = . . - , (?>2 = -— — , (p3 = ———- 
1 + ZZ I 1 + ZZ 1 + ZZ 

This defines a stereographic mapping of the plane onto the sphere 5*2, as one 
verifies by setting z = tan{r]/2) e*^, 77 G [0, tt]. 

In the 0(3) representation the topological charge 4nQ has the interpretation 
of an area of surface in 5*2, itself embedded in M"^: 



Q = ^ / dv?a A d(/?a = -^^ijk f <Pid4>j A (pk = ^t^„eijk / d^x(pid^(pjdu(pk ■ 



The result is a multiple of the area of the sphere 5*2, which in this interpretation 
explains the quantization. 
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5.3 Instantons and anomaly: non-ahelian gauge theories 

We now consider non-abelian gauge theories in four dimensions. Again gauge 
field configurations can be found that contribute to the chiral anomaly and for 
which therefore the r.h.s. of equation ( [4.50] ) does not vanish. A specially inter- 



esting example is provided by instantons, i.e. finite action solutions of euclidean 
field equations. 

To discuss this problem it is sufficient to consider pure gauge theories and the 
gauge group SU (2) since a general theorem states that for a Lie group containing 
SU (2) as a subgroup the instantons are those of the SU (2) subgroup. 

In the absence of matter fields it is convenient to use a 5*0 (3) notation. The 
gauge field is a vector that is related to the element 21^ of the Lie algebra 
used previously as gauge field by 

2t^ = -izA^-a, (5.28) 

where ai are the three Pauli matrices. The gauge action then reads: 

S{A,) = ^ I [F^.(x)f d^x, (5.29) 

{g is the gauge coupling constant) where the curvature F^j^ is also a vector: 

Ff,^ = d^A^ - d^A^ + A^ X A^ . (5.30) 

The corresponding classical field equations are 

D,F,^ = a,F,^ + A, xF,^ = 0. (5.31) 

The existence and some properties of instantons in this theory follow from con- 
siderations analogous to those presented for the CP{N — 1) model. 
We define the dual of the tensor F^^, by 

F ^1/ 2^f^iypa'^ per • (5.32) 

Then the Bogomolnyi inequality 

j d^x [f^, {x) ± F^, (x)] ' > , (5.33) 

implies: 

5(A^)>87r2|g(A^)|//, (5.34) 

with 

Q(A^) = ^ / F^, ■ F^, . (5.35) 
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The expression Q{A^) is proportional to the integral of the chiral anomaly ( [4 .431) , 
here written in SO{3) notation. We have shown that the quantity F^j, ■ F^^, is a 
pure divergence (equation (|4.44| )): 



F • F = d V 



with 



= -4 e^^p^ tr (21,9^21, + §21,21^21,) (5.36a) 
= 2e^,p, [A, ■ dpA^ + iA, ■ {Ap x A,)] . (5.366) 

The integral thus depends only on the behaviour of the gauge field at large 
distances and its values are quantized (equation (|4.50| )). Here again, as in the 



CP{N — 1) model, the bound involves a topological charge, Q{Ap). 
We now use equation ( ^.44] ) and Stokes theorem 



d'^xd^Vp = / dO fipVp , 

V JdV 

where dO is the measure on the boundary dV of the four-volume T) and ftp the 
unit vector normal to dV. We take for V a sphere of large radius R and finds 
for the topological charge Q, 

QiK) = ^ / F/^- ■ = ^/^^M : (5-37) 

The finiteness of the action implies that the classical solution must asymptotically 
become a pure gauge, i.e. with our conventions, 

^p^-\iAp-a = z{x)dpg,-\x)+0{\x\-^) |a;|^oo. (5.38) 

The element g of the SU{2) group can be parametrized in terms of Pauli matrices 

g = M4 + m ■ a , (5.39) 

where (w4, u) is a four-component real vector belonging the unit sphere S^: 

ul + u^ = 1. 

Since SU{2) is topologically equivalent to the sphere S's, the pure gauge config- 
urations on a sphere of large radius \x\ = R define a mapping from 5*3 to 5*3. 
Such mappings belong to different homotopy classes that are characterized by an 
integer called the winding number. Here we identify the homotopy group 7r3(5'3), 
which again is isomorphic to the additive group of integers Z. 
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The simplest one to one mapping corresponds to an element g(x) of the form 

g(x) = ^i^tlZi:^, r = (xl + x2)V2, (5.40) 
r 

and thus 

~ 2(x4(5,^ + e,^fcXfc)r-2, A^ = -2x,r-2. (5.41) 

r — >oo 

Note that the transformation 

g(a;) Uig(x)u5 = g(Ra;), 

where Ui and U2 are two constant SU{2) matrices, induces a SO{4) rotation of 
matrix R of the vector x^. Then 

U2a^g^(a:)U; = R^,d,g\Kx), Vig{x)d^gHx)V\ = g(Ra;)i?^,a,g^(Ra;), 
and therefore 

Introducing this relation into the definition ( |5.36a| ) of we verify that the de- 
pendence in the matrix Ui cancels in the trace and thus transforms like a 
four-vector. Since only one vector is available, and taking into account dimen- 
sional analysis we conclude 

V;, oc x^/r^ . 

For R 00 approaches a pure gauge (equation ( p.38| )), and therefore 
can be transformed into 

Vjj, ~ 3 C/xjypcr ■ (Ap X Act). 

It is sufficient to calculate Vi. We can choose p = 3, a = 4 and multiply by a 
factor six to take into account all other choices. Then 

Vi = 16eijkix4S2i + ei2iXi)(x4S3j + ej3rnXm)xk/r^ = 16xi/r^, 

and thus 

The powers of R in equation ( |5.37D cancel and since J dO = 27r^ the value of the 
topological charge is simply 

Q(A^) = 1. (5.42) 



If we compare this result with equation p.5UD we see that we have indeed found 
the minimal action solution. 
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Without explicit calculation we know already from the analysis of the index 
of the Dirac operator that the topological charge is an integer: 

Q{A^) = ^ j d^xF^,-F^,=n. (5.43) 

As in the case of the CP{N — 1) model this result has a geometric interpretation. 
In general in the parametrization ( |5.39| ), 



1 — >oo 

A few algebraic manipulations starting from 

/ R^dVL n^Vij, = \e^i,p„ I V^du„ A dup A du, 



'S3 

then yield 



Q = -T^e^upa / Updui, A dup A dua , (5.44) 



where the notation of exterior differential calculus again has been used. The area 
Tip of the sphere Sp-i in the same notation can be written 

27rP/^ _ 1 r 
^ ~ r(p/2) ~ {p- J ^Mi*^^M2 A ... A du^^ , 

when the vector Up, describes the sphere Sp-i only once. In the l.h.s. of equation 
(|5.44| ) one thus recognizes an expression proportional to the area of the sphere 
5*3. Because in general Up describes 5*3 n times when Xp describes 5*3 only once 
a factor n is generated 

The inequality (|5.35| ) then implies 

S{Ap)>8n''\n\/g\ (5.45) 

The equality, which corresponds to a local minimum of the action, is obtained 
for fields satisfying the self-duality equations 

Fp, = ±Fp, . (5.46) 

These equations, unlike the general classical field equations ( |5.31|) , are first order 
partial differential equations and therefore easier to solve. The one-instanton 
solution, which depends on an arbitrary scale parameter A, is 

2 • 2x■ 

^m = T? (^4(^*771 + Cimfca^fc) , m = l,2,3, A\ = — „ , \„ . (5.47) 
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The semi- classical vacuum. We now proceed in analogy with the analysis of 
the CP{N — 1) model. In the temporal gauge A4 = the classical minima of the 
potential correspond to gauge field components A^, i = 1,2,3, which are pure 
gauge functions of the three space variables xf. 

2im = -^i^m ■ cr = g{xi)d^g~'^{xi) . (5.48) 

The structure of the classical minima is related to the homotopy classes of map- 
pings of the group elements g into compactified (because g(x) goes to a 
constant for |x| — > 00), i.e. again of 5*3 into S'3 and thus the semi-classical vac- 
uum has a periodic structure. One verifies that the instanton solution ( |5.47|) , 
transported into the temporal gauge by a gauge transformation, connects min- 
ima with different winding numbers. Therefore, as in the case of the CP{N — 1) 
model, to project onto a ^-vacuum, one adds a term to the classical action of 
gauge theories: 

Se{A^) = S{A^) + ^Jd^x F^, ■ F^, , (5.49) 

and then integrates over all fields A^ without restriction. At least in the semi- 
classical approximation, the gauge theory depends on an additional parameter, 
the angle 9. For non-vanishing values of 9 the additional term violates CP 
conservation, and is at the origin of the strong CP violation problem, because 
if 9 does not vanish experimental bounds are consistent only with unnaturally 
small values. 



5.4 Fermions in an instanton background 

We now apply this analysis to QCD, the theory of Strong Interactions, where 
Nf Dirac fermions Q, Q, the quark fields, are coupled to non-abelian gauge 
fields A^ corresponding to the SU{3) colour group. The action can be written 
(returning here to standard SU (3) notation) 



5(A^,Q,Q) 



d^ 



X 



Nf 



The existence of abelian anomalies and instantons has several physical conse- 
quences. We mention here two of them. 

The strong CP problem. According to the analysis of Section [4.5| only con- 
figurations with a non-vanishing index of the Dirac operator contribute to the 
^-term. Then the Dirac operator has at least one vanishing eigenvalue. If one 
fermion field is massless, the determinant resulting from the fermion integration 
thus vanishes, the instantons do not contribute to the functional integral and 
the strong CP violation problem is solved. However such an hypothesis seems to 
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be inconsistent with experimental data on quark masses. Another scheme was 
based on a scalar field, the axion, which unfortunately has remained experimen- 
tally invisible. 

The solution of the U{1) problem. Experimentally it is observed that masses 
of pseudo-scalar mesons are smaller or even much smaller (in the case of pions) 
than the masses of the corresponding scalar mesons. This strongly suggests an 
approximate chiral symmetry corresponding to small quark u and d masses, 
and a more badly violated chiral symmetry corresponding to the strange quark, 
realized in a phase of spontaneous symmetry breaking. This picture is confirmed 
by its many other physical consequences. 

In a theory in which the quarks would be massless the action would have a 
chiral U{Np) x U{N-f) symmetry, in which Np is the number of fiavours. The 
spontaneous breaking of chiral symmetry to its diagonal subgroup U{N-p) leads 
to expect Np Goldstone bosons associated with all axial currents (corresponding 
to the generators of U{N) x U{N) that do not belong to the remaining U{N) 
symmetry group). If the masses of quarks are non- vanishing but small one ex- 
pects this picture to survive approximately, with instead of Goldstone bosons 
light pseudo-scalar mesons. From the preceding analysis we know that the axial 
current corresponding to the U{1) abelian subgroup has an anomaly. The WT 
identities which imply the existence of Goldstone bosons correspond to constant 
group transformations and therefore involve only the space integral of the di- 
vergence of the current. Since the anomaly is a total derivative one might have 
expected the integral to vanish. However non-abelian gauge theories have con- 
figurations that give non- vanishing values of the form (J4.5UD to the space integral 



of the anomaly (|4.43| ) . For small couplings these configurations are in the neigh- 



bourhood of instanton solutions (as discussed in section p.3|) . This indicates (but 
no satisfactory calculation of the instanton contribution has been performed) that 
for small, but non- vanishing, quark masses the U{1) axial current is far from be- 
ing conserved and therefore no corresponding light would-be Goldstone boson is 
generated. This argument resolves a long standing puzzle since experimentally 
no corresponding light pseudoscalar boson is indeed observed for Np = 2,3. 

Note that the usual derivation of WT identities involves only global chiral 
transformations, and therefore there is no need to introduce axial currents. In 
the case of massive quarks chiral symmetry is explicitly broken by soft mass 
terms and WT identities involve insertions of the operators 



Mf=mf j d^x Q/(x)75Q/ 



which are the variations of the mass terms in an infinitesimal chiral transforma- 
tion. If the contributions of A4f vanish when m/ — >^ 0, as one would normally 
expect, then a situation of approximate chiral symmetry is realized (in a sym- 
metric or spontaneously broken phase). However if we integrate over fermions 
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first, at fixed gauge fields we find (disconnected) contributions proportional to 

{Mf) = m/tr75 O + mf)~^ . 



We have shown in section 4.5 that for topologically non-trivial gauge field config- 



urations p has zero eigenmodes, which for m/ ^ give the leading contributions 

{Mf)=mfY^ I d^xip*^{x)-f5ipn{x)-^ +0{mf) 
^ J 

= (n+ — n_) -|- 0{mf). 

These contributions do not vanish for nif ^ and are responsible, after inte- 
gration over gauge fields, for a violation of chiral symmetry. 

6 Non-abelian anomaly 

We first consider the problem of conservation of a general axial current in a non- 
abelian vector gauge theory, and then the issue of obstruction to gauge invariance 
in chiral gauge theories. 

6.1 General axial current 

We now discuss the problem of the conservation of a general axial current in 
the example of an action with massless Dirac fermions in the background of 
non-abelian vector gauge fields 



S{tl;,i;;A) = - j d^x^,{x)Viip^{x). (6.1) 



In the absence of gauge fields the action S{il), ?/'; 0) has a U{N) x U{N) symmetry 
corresponding to the transformations: 

i;'= [i(l+75)U+ + |(l-75)U_]V, (6.2) 
V^' = [i(l + 75)Ul + i(l - 75)U^] , (6.3) 

where U± are N x N unitary matrices. We denote by the anti-hermitian 
generators of U{N), 

U = 1 + + 0(6*2). 

We consider gauge fields coupled to all vector currents, corresponding to diagonal 
subgroup of U{N) x U{N) corresponding to U-|- = U_, 
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We define axial currents in terms of infinitesimal space-dependent chiral trans- 
formation: 

U± = 1 ± 6'„(x)t" + 0{9^) ^6ij = 6'a(x)75t"i/', St/; = 6'a(x)i/;75t". 
The variation of the action then reads 

5S = j d'*x{J^"(a;)a^^«(a;) + ^«(a;)V5(x)757M[A^,tlV'(^)}, (6-4) 

where J^"(x) is the axial current: 

J^°(x)=^757Mt"V'- (6.5) 

Since the gauge group has a non-trivial intersection with the chiral group, the 
commutator [A^, t"] no longer vanishes 

where the fjSce-y are the totally antisymmetric structure constants of the Lie 
algebra of U{N). Thus 

5S = J d4x^«(x){-a^jf (x)-f/^«^A;^(a:)J^-^(x)}. (6.6) 

The classical current conservation equation is replaced by a gauge covariant 
conservation equation: 

D^4 = 0, (6.7) 
where we have defined the covariant divergence of the current by 

In the contribution to the anomaly the terms quadratic in the gauge fields are 
modified, compared to the expression (|4.44| ) , only by the appearance of a new 
geometric factor. Then the complete form of the anomaly is dictated by gauge 
covariance. One finds: 

DaJ^(x) = -Y^e^,p,trt-F^,F^,. (6.8) 

This is result for the most general chiral and gauge transformations. If we restrict 
both groups in such a way that the gauge group has an empty intersection with 
the chiral group the anomaly becomes proportional to trt", where are the 
generators of the chiral group G x G, and is therefore different from zero only 
for the abelian factors of G. 
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6.2 Obstruction to gauge invariance 

We now consider a non-abelian gauge field coupled to left or right-handed 
fermions for instance: 

S{i,,^-A) = - j d4xV5(x)i(l+75)P^(^), (6.9) 

(the discussion with |(1 — 75) is similar). 

The gauge theory can be consistent only if the partition function 

Z(A^) = j [diPd^] exp [-S{iP, ^; A)] (6.10) 

is gauge invariant. 

We introduce the generators t° of the gauge group in the fermion representa- 
tion and define the corresponding current as: 

j;(x)=^i(l + 75)7/.t>. (6.11) 

The invariance of Z{Kp) under an infinitesimal gauge transformation again leads 
to a covariant conservation equation for the current: 

(D^J^) =0, 

with 



= + [A 



The calculation of the quadratic contribution to the anomaly is simple: the first 
regularization adopted for the calculation in section is also suited to the 
present situation since the current-gauge field three-point function is symmetric 
in the external arguments. The group structure is reflected by a simple geometric 
factor. The global factor can be taken from the abelian calculation. It differs 
from result ( |4.21D by a factor 1/2 that comes from the projector ^(1 -|- 75). The 



general form of the term of degree three in the gauge field can also easily be 
found, but the calculation of the global factor is somewhat tedious. We argue in 
the section that it can be obtained from consistency conditions. The complete 
expression reads: 

(D^J^(x))" = -^a^e^-P- tr [t" (A^a^A, + iA^A^A,)] . (6.12) 

If the projector i(H-75) is replaced by i(l — 75) the sign of the anomaly changes. 

Unless this term vanishes identically there is an obstruction to the construction 
of the gauge theory. The first term is proportional to dap-y'- 

dc^p^ = I tr [t" (t^t^ + t^t^)] . (6.13) 
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The second term involves the product of four generators, but taking into account 
the antisymmetry of the e tensor, one product of two consecutive can be replaced 
by a commutator. Therefore the term is also proportional to dap-y. 

For a unitary representation the generators are, with our conventions, an- 
tihermitian. Therefore the coefficients d^pj are purely imaginary: 

= i tr [t- (t^t^ + t^t^)] ^ = -d^p^ . (6.14) 

For all real (the antisymmetric) or "pseudo-real" (t" = —5" '^t°'S~^) repre- 
sentations these coefficients vanish. It follows that the only non-abelian groups 
that can lead to anomalies in four dimensions are: SU{N) for iV > 3, SO{6) and 

Eq. 

6.3 Wess-Zumino consistency conditions 



In section we have calculated the part of the anomaly that is quadratic in 
the gauge field and asserted that the remaining part could be obtained from geo- 
metric arguments. The method is based on BRS transformations. The anomaly 
is the variation of a functional under an infinitesimal gauge transformation. This 
implies compatibility conditions, which here are constraints on the form of the 
anomaly. 

One convenient way to express these constraints is to express the nilpotency of 
BRS transformations. In a BRS transformation the variation of the gauge field 



takes the form: 



5brsA^(x) = D^C(x)£, (6.15) 



where C is a (fermion) ghost field and e an anticommuting constant. The corre- 
sponding variation of lnZ{A^) is: 

5BRslnZ(A^) = -/d^^ {3^ix))U^C{x)e. (6.16) 

The anomaly equation has the general form: 

(D^J^(x)) =^(A^;x). (6.17) 
In terms of A the equation ( |6.16| ), after an integration by parts, can be rewritten: 



(5BRslnZ(A^) = j d^xAiA^;x)C{x)e. (6.18) 

We then calculate the BRS variation of AC. We therefore need also the BRS 
transformation of the fermion ghost C(x): 



5BRsC(a;)=£C2(a;). 



(6.19) 
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We express that the square of the BRS operator 5brs vanishes (hke a cohomology 
operator), and thus that AC is BRS invariant 

^BRS = ^ ^BRS j d^xA{A^;x)C{x) = 0. 

This condition yields a constraint on the possible form of anomalies that deter- 
mines the term cubic in A in the r.h.s. of equation ( p.l2| ) completely. One can 
verify 

Sbrs j tr [C{x)d^, {A^dpA^ + i A^ApA^)] = . 
Explicitly, after integration by parts, the equation takes the form 

e^.pa tr J d^x {d^C^{x)A,dpA^ + d^CB^CdpA^ + dpCA^dpU^C 
+ia^C2(x) A.ApA, + ^dpC (D.CApA, + A.D^CA, + A.A^D^C) } = . 

The terms linear in A cancels automatically 

tr j d^x {dpCd^CdpA, + dpCA^dpdaC) = , 

after integrating by parts the first term and using the antisymmetry of the e 
symbol. 

In the same way the cubic terms cancel (one has to remember the anticom- 
muting properties of C) 

e;..p.tr j d4x{(a^CC + C9^C)A,ApA,+a^C([A,,C]CApA, 
+A,[Ap, C]A„ + A,Ap[A,, C])} = . 

It is only the quadratic terms that gives a relation between the quadratic and 
cubic terms in the anomaly, both contributions being proportional to 



e^i^pcT tr / d xd^Cd^CApA^ . 
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7 Lattice fermions: Ginsparg— Wilson relation 



Notation. We now return to the problem of lattice fermions discussed in section 
3 .41. For convenience we set the lattice spacing a = 1 and use for the fields the 



notation tp{x) = tpx- 

Ginsparg-Wilson relation. It had been noted, many years ago, that a poten- 
tial way to avoid the doubling problem while still retaining chiral properties in 
the continuum limit was to look for lattice Dirac operators D that, instead of 
anticommuting with 75, would satisfy the relation 

D-Ss+TsD-^ =751 (7.1) 



where 1 stands for the identity both for lattice sites and in the algebra of 7- 
matrices. More explicitly 

(D"^)^y75 + 75(D"^)^y = 75(5^y . 

More generally the r.h.s. can be replaced by any local positive operator on the 
lattice. Locality on the lattice is defined by a decrease that is at least exponen- 
tial when the point x,y are separated. The anti-commutator being local, it is 
expected that it does not affect correlation functions at large distance and that 
chiral properties are recovered in the continuum limit. Note that when D is the 
Dirac operator in a gauge background the condition (|77^) is gauge invariant. 

However, only recently have lattice Dirac operators solutions to the Ginsparg- 
Wilson relation ( fOI) been discovered, because the demands that both D and the 
anticommutator {D~^,75} should be local, seemed difficult to satisfy, specially 
in the most interesting case of gauge theories. 

Note that while the relation (|7.1|) implies some generalized form of chirality 
on the lattice, as we now show, it does not guarantee the absence of doublers, as 
examples illustrate. But the important point is that in this class solutions can 
be found without doublers. 



7.1 Chiral symmetry and index 

We first discuss the main properties of a Dirac operator satisfying relation ( [7.1| ), 
and exhibit a generalized form of chiral transformations on the lattice. 

Using the relation, quite generally true for an euclidean Dirac operator (con- 
sequence of hermiticity and refiection symmetry), 

Dt = 75D75 , (7.2) 

one can rewrite the relation ( [7.1| ), after multiplication by 75, 

D-i + (D-i)^ = l, 
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and therefore 

D + D^ = DD^ = D^D. (7.3) 
This imphes that the lattice operator D has an index, and in addition 

S = 1 - D , (7.4) 

is unitary 

SS^ = 1. (7.5) 

The eigenvalues of S lie on the unit circle. The eigenvalue one corresponds to 
the pole of the Dirac propagator. 
Note also the relations 

758 = 8+75, (758)2 = 1. (7.6) 

The matrix 758 is hermitian and ^(1 ±758) are two orthogonal projectors. If D 
is a Dirac operator in a gauge background these projectors depend on the gauge 
field. 

It is then possible to construct lattice actions that have a chiral symmetry that 
corresponds to local but non point-like transformations. In the abelian example, 

< = E(^^'"^')x,V'., (7.7) 
y 

^'^ = ^^e'^. _ (7.8) 

Indeed, 

Using the second relation ( |7.6D we expand the exponentials and reduce the equa- 
tion to 

D758 = -75D , (7.9) 

which is another form of relation ( |7.1|) . 

These transformations, however, no longer leave the integration measure over 
the fermion fields, 

]^d?/'2;d'0x , 

X 

automatically invariant. The jacobian J of the change of variables ij) \^ ij)' is 

J = dete^^^5e*^75S = 1 + 2z6'tr 75(1 - D/2) + 0(6*2), ^7_^q^ 

where trace means trace in the space of 7 matrices and in the lattice indices. 
This leaves open the possibility of generating the expected anomalies, when the 
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Dirac operator of the free theory is replaced by the covariant operator in the 
background of a gauge field, as we now show. 

Eigenvalues of Dirac operator in a gauge background. We briefly discuss the 
index of a lattice Dirac operator D satisfying the relation ( |7.1|) , in a gauge 
background. We assume that its spectrum is discrete (this is certainly true on 
a finite lattice where D is a matrix). The operator D is related by ( [7.4| ) to a 
unitary operator S whose eigenvalues have modulus one. Therefore if we denote 
by \n) its n*^ eigenvector, 

D \n) = (1 - S) \n) = (1 - e^^") \n) \n) = (1 - e"*^") \n) . 

Then using equation ( J7.2D , we infer 

D75I71) = (l-e-*^")75|7i). 



The discussion that follows then is analogous to the discussion of Section |4.5| to 
which we refer for details. We note that when the eigenvalues are not real, 9n ^ 
(mod tt), 75 |n) is an eigenvector different from |n) because the eigenvalues are 
different. Instead in the two subspaces corresponding to the eigenvalues and 
2, we can choose eigenvectors with definite chirality 

75 \n) = ±\n) . 

We call below n± the number of eigenvalues 0, and u± the number of eigenvalues 
2 with chirality ±1. 

Note that on a finite lattice Sxy is a finite matrix, and thus 

tr75(2-D) = -tr75D, 
{n\ 75(2 - D) |n) = - 5^ {n\ 75D \n) . (7.11) 



Therefore 
which implies 



In the equation all complex eigenvalues cancel because the vector |n) and 75 |n) 
are orthogonal. The sum reduces to the subspace of real eigenvalues, where the 
eigenvectors have definite chirality. In the l.h.s. only the eigenvalue contributes, 
and in the r.h.s. only the eigenvalue 2. We find 

n-i- — n_ = — (z^+ — I'-). 

The equation tells us that the difference between the number of states of different 
chirality in the zero eigenvalue sector is cancelled by the difference in the massive 
sector of eigenvalue two. 
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Remark. It is interesting to note the relation between the spectrum of D and 
75 D, which from relation ( [7. 2| ) is a hermitian matrix, 

75D = Dt75 = (75D)t, 



and thus is diagonalizable with real eigenvalues. It is easy to verify the two 
equations, the second being obtained by changing 6* in 6* + 27r, 

75D(1 - ze^^"/2 75) |n) = 2 sin(^^/2)(l -ie^'-Z^'y^) \n) , 
75D(l + ^e^^"/2 75) \n) = -2 sin(^ J2)(l + z e'^"/^ 75) \n) , 

which implies that the eigenvalues are paired ± sin(6'„/2) except for 6*^ = (mod 
tt) where \n) and 75 |?t.) are proportional. For 9n = 75D has also eigenvalue 0. 
For 6'n — TT 75D has eigenvalue ±2 depending on the chirality of \n). 
In the same way 

75(2 - D) (1 + e^^"/2 75) \n) = 2 cos(^^/2) (1 + e^^"/^ 75) \n) , 
75(2 - D)(l - e^^"/2 75) \n) = -2cos(^j2)(l - e^^"/2 75) |n) . 



Jacobian and lattice anomaly. The variation of the jacobian ( 7.1U ) can now 
be evaluated. Opposite eigenvalues of 75(2 — D) cancel. The eigenvalues 6^ = 11 
give factors one. Only 9^ — ^ gives a non-trivial contribution 



J = dete^^^s(2- 



D) 



The quantity tr75(2 — D), coefficient of the term of order 6, is a sum of terms 
that are local, gauge invariant, pseudoscalar and topological as the continuum 
anomaly (|4.43|) since 



tr75(2 - D) = ^ (n| 75(2 - D) \n) = 2(n+ 



n. 



Non-abelian generalization. We now consider the non-abelian chiral transfor- 
mations 



tl^u = + 75S)U+ + i(l - 75S)U_] ^ . 



i(l + 75)Ul + i(l-75)U 



+ 



(7.12) 
(7.13) 



where \J± are matrices belonging to some unitary group G 

u = i + e + o(e2). 
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We note that this amounts to define differently chiral components of ip and 
for -0 the definition being even gauge field dependent. 

We assume that G is a vector symmetry of the fermion action, and thus the 
Dirac operator commutes with the elements of the Lie algebra 

[D,e] = o. 

Then again the relation (|7.1J ) in the form (|7.9| ) implies the invariance of the 
fermion action 

'i/'f/ D tpu = 'i/' D -0 . 
The jacobian of an infinitesimal chiral transformation O = G+ = — G_ is 

J = 1 + tr 756(2- D) +0(e2). 



Wess-Zumino consistency conditions. To determine anomalies in the case of 
gauge fields coupling differently to fermion chiral components one can on the 
lattice also play with the nilpotency of BRS transformations, which then take 
the form 



6Cx = eC^ 



X 1 



instead of (16.151 , 16^191) . Moreover the matrix elements T^xy of the gauge covariant 
Dirac operator transform like \Jxy 

7.2 Explicit construction: Overlap fermions 

An explicit solution to the Ginsparg-Wilson relation without doublers can 
be derived from operators D-yv that share the properties of the Wilson-Dirac 
operator of equation (|3.9|), i.e. which avoid doublers at the price of breaking 
chiral symmetry explicitly. Setting 

A = 1 - Dw/M, 

where M > is a mass parameter, one takes (the idea of overlap fermions) 

S = A(AtA)"'/' ^ D = 1- A(AtA)"'/^ (7.14) 

1/2 

With this ansatz D has a zero eigenmode when A (A^ A) has the eigenvalue 
one. This can happen when A and A^ have the same eigenvector with a positive 
eigenvalue. 

Free fermions. We now verify the absence of doublers in the absence of gauge 
fields. The Fourier representation has the general form 



Dw{j>) = a{p) + i-f^f3^,{p), 



(7.15) 
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where a{p) and are real, periodic, smooth functions such that 

Pf^ip) , ?^ Pt^^ ^ 0, ?^ P^^ 

and a{p) > for all values of that correspond to doublers, i.e. such that 
(3^{p) = for \p\ ^ 0. 

In the case of the operator ( [7.15|) a short calculation shows that the determi- 
nant of D vanishes when 

1 ^ 

(M - a{p)Y + /32(p) _ M + a{p) + ^l{p) = . 

This implies P^{p) = 0, an equation that necessarily admits doubler solutions, 
and 

\M - a{p) \ = M - a{p). 

The solutions to this equation depend on the value of a{p) with respect to M 
for the doubler modes, i.e. for the values of p such that (3^{j>) = 0. If a{p) < M 
the equation is automatically satisfied and the corresponding doubler survives. 
As mentioned in the introduction, the relation ( |7. 1| ) alone does not guarantee 
the absence of doublers. If instead a{p) > M the equation implies a{p) = M, 
which is impossible. Therefore by rescaling a{p), if necessary, we can keep the 
wanted P/i = mode, while eliminating all doublers, which then correspond to 
the eigenvalue two for D, and the doubling problem is solved, at least in a free 
theory. 

In presence of a gauge field the construction can be generalized and works 
provided the plaquette action is sufficiently close to one. 

Remark. Let us stress that, if it seems that the doubling problem has been 
solved from the formal point of view, from the numerical point of view the 
calculation of the operator (A^^A)"^/^ in a gauge background represents a major 
challenge. 

8 Supersymmetric quantum mechanics and domain wall fermions 

Because the construction of lattice fermions without doublers we have just de- 
scribed is somewhat artificial, one may wonder whether there is a context in 
which they would appear more naturally. Therefore we now briefly outline 
how a similar lattice Dirac operator can be generated by embedding first four- 
dimensional space in a larger five dimensional space. This is the method of 
domain wall fermions. 

Because the general idea behind domain wall fermions has emerged first in 
another context, as a preparation, we first recall a few properties of the spectrum 
of the hamiltonian in supersymmetric quantum mechanics, a topic also related 
to the index of the Dirac operator (section |4.5|) , and stochastic dynamics or 
Fokker-Planck equation. 
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8.1 Supersymmetric quantum mechanics 

We consider a first order differential operator D that is a 2 x 2 matrix (a^ still 
are the Pauli matrices): 

D = aid^ - ia2Mx)- (8-1) 
The function A{x) is real, and thus the operator X) is anti-hermitian. 



The operator D shares several properties of the Dirac operator of section 
In particular it satisfies 

(73® + Sas = , 
and thus has an index. We introduce the operator D 

D = d^ + A{x) = -d^ + A{x), (8.2) 

and 

or 



Then 

Thus {QjQ^} are generators of the simplest form of supersymmetry, and the 
supersymmetric hamiltonian H is given by 



(8.3) 



We note that i7 is a positive operator. The eigenvectors of H have the form 
'i/'_i_(x)(l, 0) and '0_(x)(0, 1) and satisfy 

DtD|^+)=£+|^+), or DDt|V;_)=£_|V'-), £± > , (8.4) 



where 

DtD = -d^ + A^{x) - A'{x), DDt = -dl + A''{x) + A'{x). 



Moreover if x belongs to a bounded interval or A{x) oo for |a;| — > oo then the 
spectrum of H is discrete. 

Multiplying the first equation ( |8.4|) by D, we conclude that if D 7^ 0, 
and thus es^ does not vanish, it is an eigenvector of DD^ with eigenvalue £_(-, 
and conversely. Therefore except for a possible ground state with vanishing 
eigenvalue, the spectrum of H is doubly degenerate. 

This observation is consistent with the analysis of section applied to the 
operator 2). We know from that either eigenvectors are paired IV^), cral'i/') with 
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opposite eigenvalues ±i\/£, or they correspond to the eigenvalue zero and can be 
chosen with definite chirality 

2)|V') = 0, a3|V) = ±|V')- 

It is convenient to now introduce the function S{x): 

S'{x) = A{x), (8.5) 

and for simplicity discuss only the situation of operators on the entire real line. 
We assume that 

s{x)/\x\ > e>o. 

X — >±oo 

Then the function S{x) is such that e"*^*^^) is a normalizable wave function: 
J dx e~'^^^^^ < oo. In the stochastic interpretation e"^'^'^^) is the equilibrium 
distribution. 

When e~'^(^) is normalizable we know one eigenvector with vanishing eigen- 
value and chirality +1, which corresponds to the isolated ground state of D'I'D 

D|^+,0) = ^ D\^+)^0, cT3|V'+,0)-|^+,0), 

with 

On the other hand the formal solution of 'D'^\'ip-) = 0, 

is not normalizable, and therefore no eigenvector with negative chirality is found. 

We conclude that the operator D has only one eigenvector with zero eigenvalue 
corresponding to positive chirality: the index of D is one. Note that expressions 
for the index of the Dirac operator in a general background have been derived. 
In the present example they yield 

Index = I [sgn A(-|-oo) — sgn A(— oo)] , 

in agreement with the explicit calculation. 

The resolvent. For later purpose it is useful to exhibit some properties of the 
resolvent 6: 



for real values of the parameter k. Parametrizing 6 as a 2 x 2 matrix 

6 = 



Gil Gi2 
G21 G22 



61 



one obtains 

Gil = -k (D^D + k^y^ G21 = -D (D^D + A;^)"^ 
G22 = -k (DDt + Gi2 = Dt (DDt + e) . 

For k"^ real one verifies G21 = —^12- 

A number of properties tlien follow directly from the analysis of |S2 . 

When A; ^ only Gn has a pole, Gn = 0(l//c), G22 vanishes as k and 
Gi2{x,y) = —G2i{y,x) have finite limits: 

^ '^^-0^ G^2{x,y) ) 2k ^ 

Another limit of interest is the limit y —>■ x. The non-diagonal elements are 
discontinuous but the limit of interest for domain wall fermions is the average of 
the two limits 



x) = 1(1 + cr3)Gii(x, x) + i(l - a3)G22ix, x) + ia2Gi2{x, x) . 

When the function A{x) is odd, A{—x) = —A{x), in the limit a; = the matrix 
&{x, x) reduces to 

= 1(1 + ^3)^11(0, 0) + i(l - (73)^22(0, 0). 

Examples 



(i) In the example of the function S{x) = -^x^, the two components of the 



hamiltonian H become 

DD^ = -d^ + 1, D^D = -d^ -fx^ - 1. 

We recognize two shifted harmonic oscillators and the spectrum of D contains 
one eigenvalue zero, and a spectrum of opposite eigenvalues ±i^/2n, n > 1. 

(ii) Another example useful for later purpose is S{x) = \x\. Then A{x) = e(x) 
(e(x) is the sign function), and A'{x) = 26 {x). The two components of the 
hamiltonian H become 

DD^ = -d^ + l + 25(x), D^D = -d^ + 1 - 25(x). (8.6) 

Here one finds one isolated eigenvalue zero, and a continuous spectrum £ > 1. 

(iii) A less singular but similar example that can be solved analytically corre- 
sponds to A{x) = /Utanh(x), where fi is for instance a positive constant. It leads 
to the potentials 

V{x) = A\x)±A'{x)=^,'-t^^l^ 

cosh (xj 

The two operators have a continuous spectrum starting at fx^ and a discrete 
spectrum 



62 



8.2 Field theory in two dimensions 

A natural realization in quantum field theory of such a situation corresponds to 
a two-dimension model of a Dirac fermion in the background of a static soliton 
(finite energy solution of field equations). 

We consider the action "0, i/^, "0 being Dirac fermions, and ip a scalar 
boson: 

-ilj{0 + m + M^)i/j+^{d^^f + V{^) . 

We assume that V{ip) has degenerate minima, like ((/?^ — 1)^ or cos 99 and field 
equations therefore admit soliton solutions ^pix), static solitons being the instan- 
tons of the one-dimension quantum (f model. 

Let us now study the spectrum of the corresponding Dirac operator 

V = aidj: + a2dt + m + M(f{x). 

We assume for definiteness that f{x) go from —1 for x = —00 to +1 for x = +00, 
a typical example being 

f{x) = tanh(a:). 

Since time translation symmetry remains, we can introduce the (euclidean) time 
Fourier components and study 

V = aidx + iuja2 + m + M{p{x). 

The zero eigenmodes of V are also the solutions of the eigenvalue equation 

Ti \tl;) = 00 , Ti = 00 + ia2V = asd^ + ia2 (m + Mip{x)), 

which differs from equation ( |8.1| ) by an exchange between the matrices as and 
cTi. The possible zero eigenmodes of D {iv = 0) thus satisfy 

ai\il;) = e\il;), e = ±1 , 

and therefore are proportional to il^dx) solution of 

ei;'^ + {m + Mip{x))^l)^ = . 

This equation has a normalizable solution only if \m\ < \M\ and e = -fl. Then 
we find one fermion zero-mode. 

A soliton solution breaks space translation symmetry, but leads to a zero-mode 
would generate IR divergences. The zero-mode is removed by taking the position 
of the soliton as a collective coordinate and translation symmetry then is restored. 
The implications of the fermion zero- mode then require further analysis. It is 
found that it is associated with a double degeneracy of the soliton state, which 
carries 1/2 fermion number. 



dx dt 
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8. 3 Domain wall fermions 

Continuum formulation. One now considers four-dimensional space (but the 
strategy applies to all even dimensional spaces) as a surface embedded in five- 
dimensional space. We denote by the usual four coordinates, and t the coordi- 
nate in the fifth dimension. Physical space corresponds to t = 0. We then study 
the five-dimensional Dirac operator V in the background of a classical scalar field 
ip{t) that depends only on t. The fermion action reads 

S{tp, tp) = - J dt d^x ij{t, x)Vtjj{t, x), 

with 

V = + 'y5dt + Mip{Mt). 

Since translation symmetry in four-space is not broken, we can work in Fourier 
representation 

= + 75 dt + M^(Mt), 

where the parameter M is a mass large with respect to the mass of all physical 
particles. 

To find the mass spectrum corresponding to T>, it is convenient to write it 

= 7p [i\p\ + 7p75dt + 7pM(^(Mt)] , 

where 7p = p^7^/|p| and thus 7p = 1. The eigenvectors with vanishing eigen- 
value of V are also those of the operator D 

D = i'jpV + IpI = Z7p75dt + i'ypMipiMt), 

with eigenvalue |p|. 

We then note that «7p75, 7p and —75 are hermitian matrices that form a 
representation of the algebra of Pauli matrices. The operator D can then be 
compared with the operator ( |8.1J ), and M(p{Mt) corresponds to A{x). Under 
the same conditions D has an eigenvector with an isolated vanishing eigenvalue 
corresponding to an eigenvector with positive chirality. All other eigenvalues, for 
dimensional reasons are proportional to M and thus correspond to fermions of 
large masses. Moreover the eigenfunction with eigenvalue zero decays on a scale 
t = 0{1/M). Therefore for M large one is left with a fermion that has a single 
chiral component, confined on the t = surface. 

One possible interpretation of the function (p{t) is that ip{t) is a solution of 
field equations and connects two minima 99 = ±1 of the ip potential. In the 
limit of very sharp transition one is led to the hamiltonian ( ^.6|) . Note that this 
interpretation is possible only for an even dimension (i > 4; in dimension two 
like in the two-dimensional field theory, zero-modes forbid a static wall. 
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More precise results follow from the study of section |8.1| . We have noticed 
that 0(ti,t2;p), the inverse of the Dirac operator in Fourier representation, has 
a short distance singularity for t2 — > ti in the form of a discontinuity. This is 
here an artifact of treating the fifth dimension differently from the four others. In 
real space for the function t2; xi — X2) with separate points on the surface, 
xi ^ X2, the limit ti = t2 corresponds to points in five dimensions that do 
not coincide and this singularity is absent. A short analysis shows that this 
amounts in Fourier representation to take the average of the limiting values (a 
property that can easily be verified for the free propagator). Then if (fi{t) is an 
odd function one finds for ti = t2 = 

V-\p) = ^ +75) + (1 -75)/c^2(/)] , 

where di,d2 are regular functions of p^. Therefore anticommutes with 75 
and chiral symmetry is realized in the usual way. If however ifi{t) is of more 
general type one finds 

= ^ [d^ip'){l + To) + {l-l5)p'd2{p')] +ds{p'). 

where ds is regular. As a consequence 

75l)-^+I)-S5 = 2rf3b^)75, 

which is a form of Wilson-Ginsparg's relation because the r.h.s. is local. 

Domain wall fermions: lattice. We now replace four-dimensional continuum 
space by a lattice but keep the fifth dimension continuous. We replace the Dirac 
operator by the Wilson-Dirac operator ( [7.15|) to avoid doublers. In Fourier 
representation we find 

V = a(p) + if^Mli^ + 75dt + Mip{Mt). 

This has the effect of replacing by /9^(p) and shifting Mip{Mt) 1— > M(f(Mt) + 
a(p). To ensure the absence of doublers we require that for the values for which 
f^^J.{p) = ciiid p 7^ none of the solutions to the zero eigenvalue equation is 
normalizable. This is realizes if for \t\ ^ 00 (f{t) is bounded, for instance 

Mt)\<i 

and M < |a(p)|. 

The inverse Dirac operator on the surface t = takes the general form 



V-^ = i0 + 75) + (1 _ 75)52^")] + 5s{p^), 



65 



where 6i is the only function that has a pole for p = 0, and where 52,^3 are 
regular. The function ds does not vanish even if (p{t) is odd because the addition 
of a(p^) breaks the symmetry. We then always find Wilson-Ginsparg's relation 

More explicit expressions can be obtained in the limit f{t) = sgn(t) (a situation 
analogous to (|8.6|)), using the analysis of section \K^ . 

Of course to simulate domain walls on the computer one has also to discretize 
the fifth dimension. 
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APPENDICES 



Al Trace formula for periodic potentials 

We consider a hamiltonian H corresponding to a real periodic potential V{x) 
with period X: 

V{x + X) = V{x). {A\.\) 
Eigenfunctions '4>0{x) are then also eigenfunctions of the translation operator T: 

Ti;e{x) = M^ + X) = e'^^iPeix). {Al.2) 

We first restrict space to a box of size NX with periodic boundary conditions. 
This implies a quantization of the angle 9 

^iNe ^i^0^0^^ 27TP/N , 0<p<N. (A1.3) 

We call V'p,n the normalized eigenfunctions of H corresponding to the band n 
and the pseudo-momentum 9p 

NX 



and En{Op) the corresponding eigenvalues. Reality implies: 

EM = Er,i-d). (AlA) 
This leads to a decomposition of the identity operator in [0, NX] 

p,n 

We now consider an operator O that commutes with T 

[T,O] = ^ {x\0\y) = {x + X\0\y + X). 

Then 

pNX 

{q,n\0\p,m) = / dx dyi(;* ^(x) {x\0 \y) i/jp^^iy) = SpqO^niOp). 
Jo 

Its trace can be written 

pNX rX 

trO= / dx{x\0\x) = Nl dx {x\0\x) ^S^OnniOp). 
•^0 ^0 p^^ 
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We then take the infinite box hmit N ^ (yo. Then 

AT 9^ 



N ^ 2tv 

p 



d9. 







and thus we find 

r-X 1 p2Tr 



I dx {x\0\x)=y^^ [ Or,^{9)de. {A1.5) 
Jo „ 27r Jo 



We now apply this general result to the operator 

= T^e-^^ . 

Then 



which using the definition of T can be rewritten 

-X 1 r2Tt: 



[ {x + iX\e-P''\x)dx = —y[ e^^'-^""-^'^ de . {Al.7) 
Jo 27r ^ Jo 



/ [dxit)] exp = — V e'^s-pE^{e) ^ 

Jx(a/2)=x(-a/2)+ix 27r 



In the path integral formulation this leads to a representation of the form 

r.27r 

:dx(t)]exp[-5(x)] = ;^y^ / 

z{l3/2)=x{-l3/2)+eX 

where x{—P/2) varies only in [0,X], justifying the representation 
A2 Resolvent of the hamiltonian in supersymmetric QM 

The resolvent G{z) = [H + z)~^ of the hermitian operator H, 

dH = -dl + Vix), 
where —z is below the spectrum of H, satisfies the differential equation: 

(-d^ + V{x) + z)G{z; X, y) = 6{x - y) . {A2.1) 

We recall how G{z; x, y) can be expressed in terms of two independent solutions 
of the homogeneous equation 

(-d^ + +z)</Ji,2(x) = . {A2.2) 
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If we partially normalize by the wronskian W, 

W{(pi,(p2) = (p'iix)(p2ix) - ipi{x)(p2{x) = 1 , 

and moreover impose the boundary conditions 

ipi{x) for a; — s> — oo, ip2{x) — > for a; — > +00 , 
then it is easily verified that G{z; x, y) is given by 

G{z; X, y) = (pi{y)(p2ix) 9{x - y) + (pi{x)(p2iy) 0{y - x) . {A2.3) 
If the potential is an even function V{—x) = V{x), 

(P2{x) oc 

We now apply this result to the operator 

H = DD^ with z = -k^. 

The functions (fi then satisfy 

(DDt + k^) cp^ix) == [-d^ + A\x) + A'{x) + k^] cp^ix) = , 

and the equation ( |A2.3|) yields the resolvent G-{k'^; x, y). 

The corresponding solutions for the operator D^^D + k"^ follow since 

Dt (DD^ + = = (D^D + A;^) D V* = . 

Setting 

Xiix) = DVi(a^), 
we calculate the wronskian for normalization purpose 

W{xi, X2) = Xi{x)X2{x) - Xi{x)x2{x) = -k^- 
Thus the corresponding resolvent Gj^ reads 

G'+(P; X, y) = [xiiy)X2{x) e{x -y)+ Xi{x)x2{y) 0{y - x)\ . (A2.4) 

The limits x = y are 

G-{k'^]x,x) = 'pi{x)'p2{x), G+{k'^;x,x) = -^Xiix)x2ix). 
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Note that these functions satisfy third order linear differential equations. If 
the potential is even, here this implies that A{x) is odd, G±{k'^;x,x) are even 
functions. 

Finally we also need 'D^G-{k'^; x, y): 

DtG'_(fc2; y) = ^^(y)T>^^2ix)9ix -y) + ip2{y)n^M^My - x)- 

We note that D^^G_(A;^; x, y) is not continuous at x — y: 

lim B^G_{k'^;x,y) = ifi2{x)'D^ifii{x), lim DtG'_(P; x, y) = ^i(a;)DV2(a;), 

and therefore from the wronskian, 

lim DtG_(/c^;a;,y) - lim TD^G_{k'^;x,y) = 1 . 

y^x- y^x+ 

The half sum is given by 



DtG_(A;2;a;,a;) = i lim D+G_(P; x, y) + ^ lim T>^G-{k'^;x,y) 

y—>-X- ^ J/->a;+ 

= ^T>'^(fii{x)(f2ix) + ^'D^(f2ix)(fii{x) 
= \ {^1^2)' {x) + A{x)ipi{x)ip2{x). 

This function is odd when A{x) is odd. 
Finally in the limit ^ one finds 

dwe-2^("), (^2(a;) = iVe^(^) / d^^ e-^^^"), 

-oo J X 

with 

/ + 00 
du e-25(-) = 1 . 
-oo 

Moreover 

DVi(a^) = -iVe-^(^), T>^(p2{x) = ATe-^^^) . 
Therefore, as expected 

G+{k^;x,y) ~ lAr2e-s(-)-5(J/). 
fe— >o A; 

Finally 

J — oo 



and therefore 



/oo 
dt sgn(a; - t) e-^^^*) . 
-oo 



